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The  parametric  dependence  of  short  period  Lg  on  character¬ 
istics  of  the  source  ano  travel  path  is  studied  using  synthetic 
seismograms  computed  by  modal  superposition.  The  earth  models  are 
plane-layereo  ano  laterally  homogeneous  with  a  frequency-independent 
q,  and  are  meant  to  represent  the  eastern  United  states.  The  Q  and 
velocity  structure  of  the  34  km  thick  crust  are  varied  witn  the 
upper  mantle  being  represented  oy  a  halfspace. 

Synthetic  seismograms  that  qualitatively  match  many  observed 
features  of  eastern  United  states  Lg  are  computed  with  35  to  50 
mooes  in  models  that  include  some  lower  velocity  sediments  and  low  Q 
near  the  surface.  The  observations  are  represented  by  the  LRSM 
snort  period  recordings  of  the  SALMON  nuclear  explosion  in  Miss¬ 
issippi.  Three  quantitative  data  features  not  reproduced  by  the 
synthetics  are:  (1)  There  is  not  enough  energy  at  group  velocities 
less  than  2.9  km/sec;  (2)  The  synthetic  spectra  fall  off  more 
rapidly  above  2.5  Hz  than  observed  spectra  from  SALMUh;  13)  The 
amplitude  attenuation  with  range  is  much  faster  than  observed  in  the 
eastern  United  States.  While  some  of  this  may  be  due  to  not  in¬ 
cluding  the  mulcipathing  and  scattering  that  must  be  present  at 
these  high  frequencies,  the  requirement  for  frequency-dependent  Q  is 
strongly  indicated.  — . - 

Since  many  qualitative  features  of  Lg  are  reproduced  by  the 
synthetic  seismograms,  they  are  computed  for  double-couple  sources 
at  several  orientations  to  delineate  the  gross  dependence  of  Lg  on 
source  aepth.  Strike-slip  ano  45°  dip-slip  earthquakes  show  a 
similar  uniform  amplitude  decrease  with  depth.  Vertical  dip-slip 
sources  show  the  opposite  trend,  with  amplitude  increasing  with 
uepth  between  1  and  25  km  depths  in  an  eastern  United  States  model. 
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I.  introduction  and  summary 


The  objective  of  this  contract  is  to  apply  theoretical 
seismogram  methods  to  determine  tne  oependence  of  the  important 
regional  phases,  especially  Lg,  on  characteristics  of  the  source  and 
travel  path.  Our  approacn  is  to  construct  synthetic  seismograms 
that  resemole  the  data,  to  uelineate  the  parametric  dependence  of 
these  seismograms,  and  then  to  attempt  to  explain  observed  features 
of  the  data,  using  the  models  as  a  guide.  Such  an  approach  has  been 
quite  successful  in  improving  understanding  of  teleseismic  oata,  but 
has  scarcely  been  attempted  for  regional  data.  There  seem  to  De  two 
reasons  for  this.  One  is  that  the  regional  aata  characteristics 
have  not  been  studied  very  thorougnly,  out  this  deficiency  is 

rapidly  being  corrected.  The  second  is  that  the  construction  of 

realistic  synthetic  seismograms  at  regional  distances  is  relatively 
difficult  ana  convenient  techniques  have  not  been  available,  we 

have  developed  a  set  of  programs  for  superposition  of  the  normal 
modes  for  an  earth  model  which  provide  a  convenient  synthesis  tool, 
and  these  programs  are  applied  in  the  research  described  in  this 
report. 

The  main  body  of  this  report  incluoes  three  sections  ana  we 
will  summarize  the  main  results  from  each.  In  Section  II  we 

describe  our  moaal  superposition  method  for  computing  synthetic 

seismograms  and  compare  with  results  from  a  direct  wave-number 
integration  program  that  computes  the  total  solution.  Section  III 
is  concerned  with  the  synthesis  of  Lg  for  shallow  explosions.  We 
consider  crustal  models  that  represent  the  eastern  United  States  and 
study  the  effect  of  realistic  parameter  variations  on  the  character¬ 
istics  of  synthetic  Lg  compared  to  observations.  Using  one  of  the 
"best"  crustal  models  of  Section  III,  we  compute  the  theoretical 
dependence  of  Lg  on  source  depth  for  point  double-couple  repre¬ 
sentations  for  earthquakes.  The  results  are  described  in  Section  IV. 

In  this  report  we  are  only  concerned  with  synthesizing  Lg. 
All  the  earth  models  to  be  discussed  include  the  upper  mantle  as  a 
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halrspace.  Therefore,  the  cutoff  velocity  for  the  mooal  calcul¬ 
ations  is  about  4.5  km/sec,  the  shear  velocity  of  the  upper  mantle, 
ana  energy  travelling  with  greater  velocities  is  not  mcluaed  in  the 
modal  solution.  But  how  well  do  the  mooes  represent  the  total 
solution  for  apparent  velocities  less  than  4.5  km/sec?  This  can 
only  be  aetermined  (for  complex  structures)  by  comparing  with  total 
solutions. 

In  Section  II  we  compare  the  mooal  seismograms  with  seismo¬ 
grams  computed  with  the  direct  wave-number  integration  program  PROSE 
(Apsel,  1979).  The  comparison  has  some  uncertainties,  but  aoes 
indicate  that  the  modes  include  most  of  the  important  features  of 
the  seismogram  in  the  included  group  velocity  window.  More  complete 
and  unambiguous  comparisons  are  being  done  with  PkOSE  and  with  a 
direct  wave-number  integration  program  recently  developed  at 
S-Cubed,  ana  these  will  be  described  in  our  next  report. 

Most  of  the  important  results  are  included  in  our  rather 
lengthy  Section  III.  The  question  addressed  in  that  Section  can  be 
framed  as  follows:  Using  plane-layered,  laterally  homogeneous 
crustal  models  with  frequency-independent  Q,  can  the  main  observed 
characteristics  of  Lg  in  the  eastern  United  States  be  reproduced? 
The  observations  are  representea  by  the  SALMON  nuclear  explosion  in 
Mississippi.  These  observations  are  consistent  with  the  trends  ob¬ 
served  for  eastern  United  States  earthquakes,  but  are  especially 
convenient  for  synthetic  seismogram  studies  because  we  know  the 
source  function  rather  well.  We  emphasize  that  we  are  not  trying  to 
model  the  SALMON  recordings  in  detail,  but  are  using  them  to 
represent  a  broad  class  of  eastern  U.S.  Lg  observations. 

The  results  of  Section  III  can  be  quickly  summarized.  Stan¬ 
dard  earth  models  give  synthetic  Lg  seismograms  with  a  realistic 
appearance  as  long  as  they  include  some  lower  velocity,  low  Q  mate¬ 
rial  near  the  surface  (within  a  kilometer  or  so).  If  this  is  not 
included,  an  unrealistic  fundamental  model  dominates  the  seismo¬ 
gram.  we  can  probably  tnink  of  the  low  velocity,  low  Q  layers  near 
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the  surface  as  representing,  in  a  crude  way,  the  lateral  heterogene¬ 
ity  that  must  be  quite  strong  over  the  multi-hundred  kilometer  paths 
studied. 

There  are  three  basic  features  of  the  data  that  we  are  unable 
to  reproduce.  These  are: 

1.  The  syntnetics  do  not  include  the  high 
frequency  energy  arriving  with  apparent  velo¬ 
cities  slower  than  2.9  km/sec  tnat  is  a 
prominent  feature  of  the  observations. 

2.  The  synthetic  spectra  fall  off  more  rapidly  at 
high  frequencies  than  the  observed  SALMON 
spectra. 

3.  The  amplitude  attenuation  with  range  is  much 
faster  than  is  observed  in  the  eastern  United 
States  for  earthquakes  or  explosions. 

The  first  deficiency  may  be  to  some  degree  explained  by 

arguing  that  tne  late  arriving  high  frequency  energy  must  be 
associated  with  multipathing  and  scattering  by  lateral 
heterogen i ties.  This  can  never  be  reproduced  with  the  laterally 

homogeneous  models  used.  If  this  is  the  primary  reason  for  the 
first  deficiency,  we  could  focus  our  attention  on  the  portion  of  the 
wavetrain  before  2.9  km/sec  and  not  be  much  concerned  with  the  rest. 

The  second  two  deficiencies  (and,  probably,  some  of  the  first) 

are  almost  certainly  due  to  deficiencies  in  our  model  for  Q.  In 
Section  II  we  describe  some  parameter  studies  that  suggest  that  no 
frequency-indepenoent  Q  model  will  be  able  to  match  all  aspects  of 
the  data.  A  frequency-dependent  Q  is  strongly  suggested,  we  will 
oe  describing  results  from  such  models  in  our  next  report. 

The  better  crustal  models  used  in  Section  III  do  reproduce 

many  important  features  of  the  Lg  observations.  Therefore,  they 
should  be  able  to  delineate  the  gross  dependence  of  Lg  on  source 

depth.  In  Section  IV  we  compute  synthetic  seismograms  to  plot  this 

dependence.  This  is  done  for  four  sources;  an  explosion  and  point 
double-couples  at  three  basic  orientations  (strike-slip,  vertical 
dip-slip  and  45*  dip-slip).  The  explosion  source  gives  an 


unrealistic  portrayal  of  Lg  amplitude  depth-dependence,  which  is  an 
interesting  result,  even  though  explosions  are  not  practical  sources 
below  a  few  kilometers. 

The  Lg  amplitude  decreases  with  depth  at  the  same  rate  for 
strike-slip  and  45*  dip-slip  earthquake  sources.  The  decrease  is 
fairly  uniform  and  is  a  factor  of  4  to  5  between  1  ano  25 
kilometers.  The  vertical  dip-slip  source  shows  the  opposite  trend, 
with  the  amplitude  increasing  with  depth  at  about  the  same  rate. 
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II.  SEISMOGRAM  SYNTHESIS  BY  MODAL  SUPERPOSITION 


2.1  INTRODUCTION 

Our  attention  in  this  report  is  directed  towara  the  synthesis 
of  Lg  and  slower  phases.  Modal  superposition  is  an  especially  con¬ 
venient  method  for  synthesizing  this  portion  of  tne  seismogram,  we 
begin  with  a  brief  description  of  the  method.  We  will  then  compare 
synthetic  seismograms  computed  with  moaal  superposition  with  seismo¬ 
grams  computed  with  direct  wave-number  integration. 

Direct  wave-number  integration  leads  to  a  complete  solution 
which,  in  principle,  is  "exact."  However,  direct  wave-number 
integration  is  computationally  very  difficult  ano  requires  some 
numerical  compromises  that  inevitably  degrade  the  solution  to  some 
degree.  This  must  be  considered  when  evaluating  comparisons  between 
the  two  methods  for  complex  problems. 

2.2  DESCRIPTION  uF  THE  METHOD 

For  a  spherical  earth  the  normal  modes  give  the  complete 
ground  motion.  Knopoff,  Schwab  and  their  colleagues  (e.g.,  Knopoff, 
et  aU,  1973;  Panza  ano  Calcagnile,  1975;  Nakanishi,  et  al_,  1977; 
Montovani,  et  al_,  1977)  have  used  Love  wave  modal  superposition  to 
compute  synthetic  seismograms  in  numerous  studies  of  the  properties 
of  various  observed  phases.  Modes  also  give  the  complete  solution 
for  a  plane-layered  elastic  plate  over  a  rigid  halfspace.  D.  V. 
Harvey  and  C.  B.  Archambeau  of  the  University  of  Colorado  have  been 
using  this  concept  to  construct  "complete"  seismograms  for  regional 
ground  motions  from  earthquakes  (Harvey,  1980). 

We  have  developed  a  program  based  on  Knopoff 's  method  (Schwab 
and  Knopoff,  1972)  that  allows  rapid  and  inexpensive  computation  of 
all  Rayleigh  ano  Love  wave  modes  up  to  the  highest  frequencies  of 
interest  for  seismological  problems.  The  root  finding  program 
improves  on  the  basic  Schwab  and  Knopoff  algorithm  by  taking 
advantage  of  the  frequency  independence  of  many  of  the  computed 
parameters  as  the  mode  search  is  conducted  at  a  fixed  phase  velocity. 
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The  published  Schwab  and  Knopoff  (1972)  algorithm  is  for 
calculation  of  the  Love  and  Rayleigh  secular  functions  only.  For 
seismogram  synthesis,  excitation  parameters  and  modal  eigenfunctions 
are  required.  The  Laplace's  expansion  technique  used  by  Schwab  and 
Knopoff  has  been  extended  to  compute  these  quantities  using  the 
interval  elements  of  the  root-finding  program.  The  computed  eigen¬ 
functions  seem  to  be  relatively  insensitive  to  the  instaDilities 
that  compromise  the  accuracy  of  high  frequency  calculations  with 
previously  available  surface  wave  programs. 

2.3  COMPARISON  OF  MOUAL  SUPERPOSITION  WITH  THE  COMPLETE  SOLUTION 

FROM  PROSE 

For  a  standard  earth  model,  the  normal  modes  represent  only  a 
portion  of  tne  solution.  But  the  excluoeo  part  of  the  solution  may 
not  be  very  important  for  many  problems  of  interest.  Unfortunately, 
it  is  not  easy  to  know  what  the  modal  solution  is  missing  without 
doing  the  full  calculation  to  see.  In  this  section,  we  will  compare 
some  modal  solutions  with  total  solutions  obtained  with  direct 
wave-number  integration. 

The  layered  earth  model  cnosen  for  our  comparison  calculations 
is  listed  in  Table  1.  The  first  50  Rayleigh  wave  modes  were  com¬ 
puted  for  this  structure  and  the  phase  ana  group  velocity  dispersion 
are  plotted  in  Figure  1.  The  cutoff  velocity  is  4.52  km/sec,  the 
shear  velocity  of  the  halfspace.  The  modal  solution  can  include  no 
energy  traveling  faster  than  tnis  velocity. 

The  comparison  calculations  were  aone  oy  ur.  John  urcutt, 
bcripps  Institute  of  Oceanography,  with  the  PROSE  program  (Apsel, 
1979;  Apsel,  et  aj_. ,  1930).  This  program  is  much  more  expensive 
than  the  modal  superposition  program  and  the  cost  scales  up  rapidly 
with  frequency  ana  source-receiver  range.  The  problem  parameters 
were  chosen  to  minimize  the  expense. 

Synthetic  seismograms  were  computed  at  four  ranges  (50,  100, 
150,  200  km)  for  a  point  double-couple  representing  strike-slip 
faulting  at  a  depth  of  1  km.  The  azimuth  was  45  degrees  from  the 
strike  and  the  source  time  function  was  a  step  with  a  moment  of 
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TABLE  1 


t. 

EARTH  MODEL  FOR 

COMPARING  MODE  AND  DIRECT  WAVE-NUMBER 

INTEGRATED  SEISMOGRAMS 

Depth 

(km) 

Thickness 

(km) 

a 

(km/sec) 

s 

(km/sec) 

p 

(gm/sec3) 

t; 

1.5 

1.5 

3.7 

2.16 

2.1 

35 

8.0 

6.5 

6.1 

3.3 

2.85 

250 

. 

34.0 

26.0 

6.6 

3.59 

3.05 

1800 

c 

t. 

00 

00 

8.1 

4.52 

3.35 

2000 

• 
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10  dyne-cm.  The  WwShN  short  period  instrument  response  was 
applied  and  the  calculations  were  done  with  a  Nyquist  frequency  of  5 
Hertz  after  tapering  the  spectrum  to  zero  between  0.5  and  1.0  Hertz 
with  a  cosine-squarea  filter.  (Only  thirteen  modes  contribute  to 
frequencies  less  than  1  Hertz.;  waking  appropriate  adjustments  for 
the  fact  that  the  spectrum  is  tapered  after  applying  the  seismometer 
response,  the  tapering  filter  could  be  considered  part  of  the  source 
spectrum.  It  has  a  strong  influence  on  the  solution. 

The  seismograms  are  compared  in  Figure  2  in  two  ways,  on  the 
basis  of  their  waveforms  and  on  the  basis  of  their  peak-to-peak 
amplitudes.  The  waveform  comparisons  show  that  the  modes  include 
most  of  the  important  features,  especially  for  group  arrival  times 
associated  with  velocities  of  3.5  km/sec  and  slower.  Energy 
arriving  with  group  velocities  much  faster  than  3.5  km/sec  is  very 
small  for  this  example. 

The  amplitude  comparison  is  not  very  satisfactory.  We  suspect 
some  minor  incompatibility  in  our  specification  of  the  filters  used 
to  construct  the  seismograms,  and  are  attempting  to  resolve  it.  We 
have  computed  seismograms  with  both  methods  in  a  halfspace  where  the 
exact  Rayleigh  wave  solution  is  easily  computed.  The  methods  agree 
within  a  few  percent  when  the  true  ground  motion  (no  filters)  is 
compared. 

An  irresolvable  difference  between  the  PROSE  ana  mode  calcu¬ 
lations  is  in  the  specification  of  attenuation.  In  PROSE,  attenu¬ 
ation  is  introduced  by  computing  with  complex  velocities.  The  route 
could  be  followed  for  our  modal  calculations,  but  only  with  a  large 
increase  in  computational  effort.  we  chose  to  include  Q  by 
computing  an  effective  attenuation  factor  using  the  theory  intro¬ 
duced  Dy  Anderson  and  Archambeau  (1964).  This  first-order  pertur¬ 
bation  theory  is  usually  used  for  surface  waves.  The  two  methods 
should  give  nearly  the  same  results  for  high  Q,  but  our  structure 
(Table  1)  includes  low  Q  in  the  near  surface  layers.  Our  impression 
is  that  most  of  the  waveform  and  amplitude  differences  oetween  the 
PROSE  and  modal  seismograms  in  Figure  2  is  due  to  differences  in  the 
way  the  Q  is  introduced. 


16.65 


PROSE 


TIME  (s) 


Figure  2.  Synthetic  seismograms  computed  with  modal  superposition 
are  compared  with  those  computed  with  PROSE.  Zero  on 
the  time  scale  of  the  plots  Is  reduced  by  R/8.1.  The 
maximum  peak-to-peak  amplitude  in  microns  Is  shown  at 
the  left  and  several  Interesting  group  velocities  are 
marked  in  each  set. 


III.  SYNTHESIS  OF  lg  FOR  SHALLOW  EXPLOSIONS 


3.1  INTRODUCTION 

The  Lg  phase  is  a  high  frequency  wavetrain  which  arrives  with 
a  group  velocity  near  3.5  km/sec.  It  is  observed  on  all  three 
components  for  both  explosions  and  earthquakes.  In  this  section,  we 
address  the  following  question;  using  modal  superposition  for  plane¬ 
layered  earth  models,  can  realistic  Lg  phases  be  synthesized  for 
shallow  explosions? 

To  auaress  this  question  we  consider  only  models  that  are  ap¬ 
propriate  for  the  eastern  Uniteo  States  (EUS).  The  gross  character¬ 
istics  of  the  crust  in  this  region  are  fairly  well  unoerstood,  but 
details  vary  from  place-to-place.  We  compute  seismograms  for  a  num¬ 
ber  of  different  models,  but  they  are  all  meant  to  be  "appropriate 
for  the  EUS." 

One  reason  for  focussing  on  EuS  crustal  models  is  that  this 
region  is  representative  of  stable  platform  areas  that  also  occur  in 
Eurasia.  Also,  there  was  a  nuclear  explosion  experiment,  SALMON 
(5.3  KT)  near  the  southern  edge  of  the  region.  This  event,  which 
was  at  a  depth  of  830  meters  in  Mississippi,  was  recorded  at  many 
stations  in  the  EUS,  The  Lg  phase  on  these  recordings,  though  not 
as  clear  as  for  many  earthquakes,  serves  as  a  convenient  and  useful 
standard  for  comparison. 

before  continuing,  we  should  say  what  it  means  to  synthesize 
Lg  seismograms  that  "match"  the  data.  We  will  oe  using  the  LRSM 
short  perioo  recordings  of  SALMON  (see  Figure  20  in  Section  3.4)  to 
represent  the  data,  since  it  is  the  only  large  explosion  in  the 
eastern  U.S.  and  does  have  a  source  function  that  is  rather  well 
known.  However,  the  Lg  phase  is  not  very  prominent  on  the  LRSM 
short  period  seismograms  and  the  onset  is  not  easy  to  time. 
Qualitatively,  the  Lg  is  a  high  frequency,  extended  wavetrain  that 
has  a  group  arrival  time  near  3.5  km/sec.  As  we  study  the  phase  in 
more  detail,  quantitative  features  that  become  important  are  the 
approximate  duration,  the  spectral  asymptotes  and  the  amplitude. 
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We  find  that  syntnetic  seismograms  which  include  Lg  phases 
that  are  qualitatively  similar  to  the  observations  can  be  computed 
witn  earth  models  that  include  sufficient  detail.  However,  simple 
moaels  deauceo  from  observations  of  the  fundamental  mode  surface 
waves  do  not  give  realistic  seismograms.  It  is  necessary  to  have 
some  low  velocity,  low  Q  material  at  the  surface.  Gradients  in  the 
mid-crust  seem  to  be  required  if  the  arrival  is  to  be  emergent,  as 
it  often  is.  On  the  other  hand,  sharp  Lg  arrivals  can  be 

synthesized  with  a  tnick,  uniform  layer  at  the  base  of  the  crust. 

We  also  conclude  that  the  slow  (<2.9  km/sec)  high  frequency 

arrivals  seen  on  the  observed  records  are  probaoly  impossible  to 

compute  with  plane-layered  models.  Some  scattering  and/or  multi- 
pathing  is  probably  required  to  reproduce  this  portion  of  the  data. 

To  develop  the  conclusions  listed  in  the  last  two  paragrapns, 
we  proceed  in  the  following  way.  We  oegin  with  synthetic  seismo¬ 

grams  for  very  simple  models.  They  turn  out  to  be  unlike  the  data, 
ay  trial-ana-error  we  add  required  layers  of  complexity  until  the 
synthetics  have  a  realistic  appearance,  witn  the  SALMON  observations 
serving  as  a  ^qual itati ve)  standard  for  comparison.  Parametric 
studies  are  then  conducted  to  identify  the  important  controlling 
parameters,  we  emphasize  that  we  are  not  trying  to  model  records  in 
detail,  but  to  reproduce  gross  features  of  observed  Lg. 

An  important  tool  for  understanding  the  controlling  features 
of  the  models  is  introduced  in  Section  3.2.  This  is  a  graphical 
display  of  the  kinetic  energy  distribution  for  particular 
frequencies  as  a  function  of  layer  and  mode.  With  tnis  display  we 
can  see  where  the  energy  is  trapped  and  which  modes  are 
contributing.  The  display  is  presented  for  a  model  that  has  most  of 
the  required  features. 

Using  the  knowledge  gained  from  the  models  described  in 
Section  3.2,  we  constructed  our  preferred  EUS  crustal  model,  called 
SI.  This  model  is  introduced  at  the  beginning  of  section  3.3. 
Synthetic  seismograms  are  compared  to  the  SALMON  data  in  that 
section  in  both  the  time  and  frequency  domain.  Our  conclusions 


about  this  comparison  are  summarized  at  the  eno  of  Section  3.3. 
basically,  the  synthetics  are  qual itatitively  like  the  data  in  the 
group  velocity  winaow  from  3.6  to  2.9  km/sec  and  in  the  frequency 
Dana  from  0.5  to  2.5  Hz  (where  the  data  have  large  signal/noise). 

The  synthetics  for  the  model  SI  fail  to  account  for  late 
arriving  (slower  than  2.9  km/sec)  high  frequency  energy  that  is 
prominent  on  the  observea  records.  Also,  frequencies  aoove  2.5  Hz 
are  poorly  modeleo,  especially  at  the  close  (a  «  24U  km)  stations. 
Also,  in  Section  3.4  we  examine  the  amplitude  attenuation  for  the 
model  ana  find  that  it  is  much  too  large  compared  to  eastern  United 
States  ooservations. 

In  Section  3.5  we  describe  some  parametric  studies,  primarily 
of  the  Q  model,  carried  out  to  determine  whether  the  model 
deficiencies  can  be  remedied.  Our  parametric  studies  are  far  from 
complete,  but  strongly  support  the  conclusion  that  plane-layered, 
constant  Q,  laterally  homogeneous  models  are  not  adequate  for 
accounting  for  those  model  deficiencies  described  in  the  last 
paragraph.  Our  conclusions  regarding  the  synthesis  of  Lg  for 
shallow  explosions  are  summarized  in  more  detail  in  Section  3.6. 

3.2  Lg  IN  SIMPLE  CRUSTAL  MODELS 
3.2.1  Introduction 

Tnis  is  basically  an  introuuctory  section  in  which  we  describe 
some  of  our  first  synthesis  efforts.  The  results  of  these 
calculations  lead  us  to  the  more  realistic  anc  useful  parameter 
studies  described  in  subsequent  sections.  However,  we  describe 
these  initial  modeling  efforts  as  a  reminder  of  what  must  be 
included  to  synthesize  realistic  Lg  seismograms. 

The  important  points  we  wish  to  make  in  tnis  section  may  be 
briefly  summarized.  Our  first  synthetics  were  done  for  a  simple 
three-layer  model  (Table  2)  which  is  roughly  consistent  with  the 
gross  charcteristics  of  the  crust  in  the  eastern  United  States  as 
inferred  from  fundamental  mode  dispersion  and  refraction  data.  The 
Q  model  is  based  on  that  of  Herrmann  and  Mitchell  (1975). 


All  of  the  models  discussed  in  this  report  are  Dased  on  the 
model  in  Table  1.  We  ado  details  and  make  changes  when  they  seem 
required,  but  attempt  to  hold  the  gross  features  constant.  In 
particular,  the  crustal  thickness  and  the  velocities  in  the  vicinity 
of  the  crust-mantle  boundary  are  not  changed. 

The  initial  seismograms  are  not  realistic  because  they  are 
dominated  by  a  pulse-like  fundamental  mode.  Thus,  our  first 
correction  is  to  add  lower  velocity  sedimentary  layers  at  the 
surface  to  disperse  this  mode,  we  find  that  the  Rayleigh  wave  is 
still  too  large.  Hence,  we  are  led  to  include  a  very  low  Q  in  the 
top  500  to  1000  meters  to  damp  the  unrealistic  contribution  from 
this  mode.  These  two  ideas,  the  requirement  for  low  velocity  and 
low  Q  in  the  very  near  surface  layers,  are  held  rather  consistently 
throughout  this  report. 

There  are  some  other  interesting  ideas  introduced  in  this 
section.  It  is  tempting  to  think  that  the  main  trends  can  be 
deduced  by  computing  synthetic  seismograms  with  relatively  few 
mooes,  we  show  that  this  is  not  the  case,  but  that  twenty  or  more 
modes  are  requirea. 

We  also  introduce  the  comparison  of  syntnetic  Lg  seismograms 
with  observations  of  the  SALMON  explosion.  The  LRSM  short  period 
recordings  of  tnis  event  are  plotted  in  Figure  20  in  Section  3.4. 
In  subsequent  sections,  we  will  make  extensive  use  of  these  data  as 
a  standard  to  judge  the  realism  of  our  synthetics. 

Near  the  end  of  this  section  (sub-section  3.2.7)  we  introduce 
an  important  tool  for  understanding  the  controlling  features  of  the 
crustal  models.  This  is  a  graphical  display  of  the  kinetic  energy 
distribution  for  particular  frequencies  as  a  function  of  layer  ana 
mode.  With  this  display  we  can  see  where  the  energy  is  trapped  and 
which  moues  are  contributing. 

Finally,  at  the  end  of  tnis  section  we  list  the  conclusions 
drawn  from  this  initial  investigation  of  modal  Lg  synthetics.  The 
models  introduced  in  our  parametric  studies  in  later  portions  of  tne 
report  are  motivated  in  large  part  by  these  conclusions. 
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TABLE  2 

SIMPLE  EASTERN  UNITED  STATES  CRUSTAL  MUDEL 


Uepth 

(km) 

Thickness 

(km) 

a 

(km/ sec) 

8 

(km/sec) 

p 

(gm/cc) 

Q 

7 

7 

6.1 

3.3 

2.85 

250 

17 

lu 

6.6 

3.59 

3.05 

250 

34 

17 

6.6 

3.59 

3.05 

00 

00 

8.1 

4.52 

3.35 

J 
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Fifty  mooes  were  computed  for  the  model  in  Table  2.  The 
dispersion  curves  are  plotted  in  Figure  3.  Synthetic  seismograms 
were  computed  at  several  ranges  for  an  explosion  (reduced 
displacement  potential)  source  at  several  deptris.  The  seismograms 
at  a  range  of  400  km  arid  1,  3,  5,  10  km  depths  are  shown  in  Figure 
4.  The  WWSSN  short  period  instrument  response  is  included. 

The  seismograms  appear  reasonable  at  the  deeper  depths. 
However,  at  shallow  depth  they  are  dominated  by  what  is  essentially 
a  Kayleigh  wave  in  a  halfspace  (the  surface  layer;.  At  the  1  km 
depth  this  phase  travels  at  about  3.05  km/sec  and  has  a  period  near 
1  Hz.  That  is  on  the  flat  part  of  the  fundamental  mode  group 
velocity  curve.  It  is  clearly  necessary  to  disperse  this  wave  if  we 
are  to  compute  realistic  seismograms  for  shallow  sources.  This  can 
be  done  by  adding  some  lower  velocity  material  near  the  surface. 

3.2.3  Addition  of  Surface  Sediments  to  the  Model 

The  same  model  with  some  sedimentary  layers  added  at  the 
surface  is  listed  in  Table  3.  These  are  like  the  sedimentary  layers 
inferred  for  the  Cincinnati  Arch  region  by  Herrmann  (1969)  using 
surface  wave  observations.  The  dispersion  curves  for  50  modes  are 
shown  in  Figure  5.  The  group  velocity  curves  are  not  much  different 
from  those  in  Figure  3  except  for  the  appearance  of  occasional  modes 
that  crop  below  3  km/sec. 

oei smograms  were  computed  for  the  same  range  and  depths  as  for 
the  simpler  structure  and  are  shown  in  Figure  6.  A  fairly  realistic 
Lg  phase  now  appears.  However,  the  fundamental  mode,  though  no 
longer  dominating,  shows  up  for  the  shallow  sources  as  a  damped 
sinusoid  that  is  quite  unlike  real  data.  To  eliminate  this 
unrealistic  feature,  it  is  necessary  to  nave  very  low  Q  in  the  near 
surface  region. 

3.2.4  Addition  of  Low  Q  Near  tne  Surface 

In  Table  4  we  list  a  new  model  which  is  identical  to  the  model 
of  Table  3  below  2.8  km.  8ut  we  have  now  added  a  larger  velocity 


TABLE  3 


EASTERN  UNITED  STATES  CRUSTAL  MODEL 
WITH  SURFACE  SEDIMENTS 


Depth 

(km) 

Thickness 

(km) 

a 

(km/ sec) 

B 

(km/sec) 

p 

igm/cc) 

Q 

0.6 

u.6 

3.7 

2.16 

2.50 

200 

4.1 

3.5 

5.6 

3.14 

2.67 

200 

7.0 

2.9 

6.1 

3.3 

2.85 

250 

17.0 

10.0 

6.6 

3.59 

3.05 

250 

34.0 

17.0 

6.6 

3.59 

3.05 

2000 

00 

00 

8.1 

4.52 

3.35 

2000 
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The  phase  and  group  velocity  dispersion  Is  plotted  for  the  first  50  modes  of  the 
four-layer  eastern  United  States  crustal  model  In  Table  2. 
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Figure  4.  Synthetic  seismograms  are  shown  at  several  depths  for  the  simple 
eastern  United  States  crustal  model  of  Table  2,  The  source  is  a 
reduced  displacement  potential  with  a  corner  frequency  near  1  Hz 
and  the  range  Is  400  km.  Fifty  modes  were  summed  and  the  WWSSN 
short-period  Instrument  response  Is  included.  The  peak  amplitude 
In  microns  at  1  Hz  Is  Indicated  t  the  left  or  each  record.  The 
time  of  a  3.5  km/sec  arrival  Is  marked.  Zero  on  the  time  scale 
Is  the  origin  time  +  30  seconds. 


simple  eastern  United  States  crustal  model  with  sedimentary  layers  added  at  the 
surface.  (Table  3) 
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TABLE  4 


EUS  CRUSTAL  MODEL  S3 


Depth 

Thickness 

a 

B 
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gradient  in  the  top  2.3  km  and  a  very  low  Q  in  the  top  600  meters. 
The  phase  ana  group  velocity  curves  for  the  first  35  modes  are  shown 
in  Figure  7. 

In  Figure  8  we  compare  a  20  mode  synthetic  seismogram  to 
recordings  of  the  SALMON  event.  The  range  for  the  synthetic  is  244 
km  compared  to  244.6  km  (EUAL)  and  243.5  km  (OELA)  for  the 
observations.  The  amplitude  comparison  is  not  meaningful.  The 
explosion  depth  for  the  synthetic  is  1  km  compared  to  tne  0.83  km 
depth  of  SmLMON.  We  see  that  the  synthetic  is  too  short  in  duration 
and  is  too  long  period.  The  highly  dispersed,  late  arriving  long 
period  waves  are  fundamental  mode  and  simply  are  not  seen  in 
ooservations,  though  these  ooserved  records  are  cut  off  too  soon  to 
demonstrate  that  point. 


To  be  able  to  change  the  model  to  more  closely  resemble  the 
observations,  we  need  to  understand  the  effects  of  the  controlling 
parameters.  In  Figure  9,  we  compare  the  synthetic  from  Figure  8  to 
synthetics  from  two  modified  versions  of  the  earth  model.  The  model 
for  the  three  synthetics  differs  only  in  the  depth  to  which  Q  is 

P 

assumea  to  be  20.  In  the  top  recora  this  depth  is  600  meters  (Table 
2),  in  the  middle  it  is  1.6  km  and  at  the  bottom  it  is  2.8  km.  The 
amplitudes  of  the  Lg  do  scale  down  with  the  lower  Q  ,  but  not 

P 

nearly  so  much  as  those  of  the  fundamental  mode. 


3.2.5  The  Effect  of  Higher  Modes 


In  Figure  10,  we  compare  5  mode  and  20  mode  seismograms  for 
two  cases.  The  importance  of  the  higher  modes  is  clearly 
displayed.  The  amplitude  of  the  high  frequency  waves  increases  by 
about  a  factor  of  four  when  the  higher  modes  are  added.  The 
prominence  of  the  funaamental  mode  is  correspondingly  reduced. 


3.2.6  A  Comparison  with  SALMON  Observations 


How  well  do  the  synthetic  seismograms  for  this  model  compare 
with  data  at  larger  ranges?  In  Figure  11  we  compare  two  SALMON 
recordings  with  two  synthetic  seismograms  at  nearly  the  same  range. 
The  signal /noise  ratio  is  not  very  high  for  the  FOTX  recording,  so 
it  should  not  be  taken  too  seriously  as  a  standard  of  comparison. 
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Figure  9.  The  synthetic  seismogram  for  R  =  244  km  from  Figure  8  is  compared  t< 
synthetics  from  two  modified  versions  of  the  same  earth  model  The 
seismograms  are  computed  with  20  modes. 


5  modes 


Figure  10.  The  effect  of  higher  modes  Is  shown  for  two  structures.  The  model  for  the 
top  pair  has  Q  =  20  to  1.6  km  and  that  for  the  bottom  pair  has  Q  =  20  to 
2.8  km.  The  20  mode  seismograms  are  repeated  from  Figure  9. 
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Figure  11.  Comparison  of  observed  SALMON  recordings  to  synthetic  seismograms  for  two  Q 

models.  Note  the  range  difference.  The  times  corresponding  to  arrivals  at  3.5 
and  2.9  km/sec  are  marked  on  the  seismograms. 


The  depth  of  the  Q  *  20  layer  is  the  only  feature  distinguishing  the 
synthetics.  The  model  with  the  smaller  low  Q  zone  seems  preferable 
because  of  its  longer  duration.  It  would  appear  to  be  qualitatively 

similar  to  the  VUlU  observation  if  it  were  shifted  to  later  time. 

However,  even  then  there  is  too  little  high  frequency  energy  in  the 
synthetic. 

3.2.7  Kinetic  Energy  Distribution  Plots 

The  group  velocity  curves  in  Figure  7  are  some  help  in 

understanding  what  features  of  the  model  are  controlling  the 

synthetic  seismograms.  But  what  is  really  needed  is  an  indication 
of  where  the  important  energy  is  being  trapped.  An  excellent  way  to 
obtain  this  is  to  look  at  the  kinetic  energy  for  particular 

frequencies  as  a  function  of  layer  and  mode  number. 

A  tabular  portrayal  of  the  kinetic  energy  characteristics  of 
the  moael  is  presented  in  Table  5.  In  this  table  we  list  the 
quantities  K^,  and  r.  for  three  frequencies,  I,  2,  and  3 

Hz.  The  b,  ana  depth  to  the  bottom  of  each  layer  are  also 

shown  with  each  quantity.  The  K..,  E--  and  F  are  computed 

*  J  '  J  ' 

the  following  way: 

1.  A  target  frequency  (f)  is  chosen,  in 

this  case  1,  2,  ana  3  Hz. 

2.  The  root  finaing  algorithm  varies  f 
for  fixed  phase  velocity  (c).  The 
closest  f  to  the  target  frequency  is 
found.  For  each  mode  this  frequency 
and  the  associatea  phase  and  group 
(U)  velocities  are  printed.  These 
are  the  values  plotted  in  Figure  7. 

3.  For  each  selected  f  we  compute 
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where  1^  denotes  tne  kinetic  energy  in  tne  ith  layer 
for  the  jth  mode. 

4.  Also  computed  is 
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where  Ap(j)  is  the  modal  amplification  and 
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J  u.  Q. 
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where  Q.  '  is  the  effective  Q„  for  the  jth  mode. 
J  B 

5.  Finally,  we  compute 


m  a 

a  v'4  r  i 

Tj  a  2-j  E. .  and  r.  =  -  x  100 

1  7"?  ij  i  max  . 

J*1  1«l,n  1 

Tnus,  for  each  selected  frequency  we  print  three  tables.  The 
first  displays  the  distribution  of  kinetic  energy  for  each  mode. 
For  the  second  table  the  energies  are  weighted  by  the  factors  useo 
to  compute  synthetic  seismograms  (excluding  the  source  depth 
dependent  eigenfunctions).  These  weighted  energies  are  normalized 
to  allow  comparison  from  mode-to-mode.  Finally,  the  weighted  energy 
in  each  layer  is  summed. 

A  graphical  display  makes  it  easier  to  absorb  the  information 
contained  in  tables  like  those  in  Table  5.  In  Figure  12  we  show 
plots  of  the  K.j  and  E . ^  distribution  for  four  cases.  The  mode 
number  versus  depth  tables  are  displayed  by  shading  the  appropriate 
positions  in  a  two-dimensional  grid,  with  the  degree  of  shading 
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TABLE  5.  KINETIC  ENERGY  DISTRIBUTION  FOR  HOOEL  S3. 


m 

r 

A 

a 

. 

: 

ft 

c 

7 

• 

9 

i 

awi 

i.’: 

.a: 

•  * 

•  • 

.  . 

.? 

.C 

.c 

2 

a. : 

!  .c 

*.i 

:i.? 

Z2.2 

9.1 

a 
.  •« 

A 
•  » 

.? 

j 

i..r 

*•*. 

‘-.5 

•  a  » 

i:.» 

'.2.2 

**.» 

.C 

• 

•  - 

* 

*  •  -  • 

r.t» 

1  -  •  4 

•  a  a 
»  •  .  • 

:.? 

.: 

* 

a.i 

!  •  ‘ 

•  c 

.ft 

.* 

1.2 

iC.I 

.0 

* 

.9* 

*.6! 

!.«! 

r.» 

a: 

•  • 

■  •  . 

1.2 

2.2 

:ai 

55.1 

•  1 

f 

r.«? 

:.r 

i.o 

3.’ 

1.4 

a: 

.as 

45. « 

•  1 

i 

a;: 

ft.* 

■*#4 

a: 

4.: 

:a* 

«... 

.2 

« 

:.3t 

2.2* 

!.* 

!.t 

ft.* 

22. 1 

•  M 

.3 

.«« 

z.n 

i.i; 

J.* 

2.2 

as 

S.’ 

:*..* 

SC. 5 

•  9 

a 

TP  » 

?.? 

i.» 

•  9 

.  . 

«•: 

!•.! 

52.2 

ai 

:: 

.»♦ 

».i: 

!.• 

a; 

•  9 

2.» 

t.c 

:m 

SI.* 

i.e 

12 

•  «c? 

’.c: 

x..« 

•  •  ■ 

a: 

9*9 

2ft.« 

•9.1 

7.9 

!• 

•  •  A  4 

•••£ 

*.3s 

a? 

i.' 

i.i 

».i 

2.o 

6.4 

22.? 

*1.2 

i*.i 

* 

f 

o 

J 

. 

■ 

*• 

5 

4 

7 

4 

• 

l 

as* 

l.*2 

1.21 

64.7 

31.9 

1.7 

.  - 

.  « 

A 

.r 

.e 

a’" 

2.5. 

an 

0.5 

1?.7 

**.7 

••1 

•  z 

.2 

v  7 

.  z 

.3 

t 

a.i 

:.t« 

!.ll 

5i.i  • 

I.I 

12.9 

.1 

.« 

•  C 

.e 

•  3 

A 

•  •  m  • 

?.*: 

2.2  2 

a.: 

!.• 

-•.9 

7.1 

•  1 

.? 

•  4. 

.3 

r 

2.22 

2.4. 

-.2 

2.9 

•  A  a 

m  •  •  • 

»«  .* 

12.7 

.2 

.5 

.  <J 

t 

..Ji 

1.22 

!.« 

1.2 

m.r 

5.2 

3.2 

72.: 

2.3 

•  3 

.3 

7 

2.21 

:.!! 

?.=♦ 

1.7 

1.3 

5.1 

5.1 

•  *  .c 

*7.: 

25.3 

•  r 

.2 

a 

•  .  a — 

T.*» 

2.3* 

... 

1.1 

3.- 

2.4 

2.3 

53.3 

2C.1 

•  s 

4 

a** 

?.*: 

3. *4 

'  .1 

... 

•  2 

.2 

1.3 

14.4 

if. 3 

.3 

1* 

..»! 

*.»i 

5.54 

•  . 

•  2  . 

•  2 

.2 

.» 

1.4 

32.4 

44.3 

•  3 

11 

2.42 

2.52 

• 

.  - 

•  1 

»• 

.4 

.»■ 

2.2 

33 .2 

5*. • 

•  V 

.  * 

.  — 

2.4* 

!.*i 

.2 

.1 

•  9 

•  3 

1.1 

23.2 

44.4 

.0 

12 

2.31 

3.-4 

3... 

•  2 

.2 

.9 

1.2 

1  .6 

0.9 

32.3 

54.4 

.3 

l» 

l.*« 

2. f. 

S.2J 

.3 

•  3 

1.2 

l.» 

2.2 

2.4 

33.3 

44.1 

.1 

1! 

av 

Z.li 

3.3. 

••  - 

.9 

1.2 

1.3 

2.2 

ft.» 

12.2 

55.2 

.1 

14 

2..  2 

’.'* 

2.11 

.5 

,4 

1.4 

2.1 

:.* 

!.! 

32.4 

54.1 

.1 

12 

2.21 

a-s 

Z.2i 

.5 

•  4 

1.4 

2.1 

2.4 

5.7 

21.3 

5*. 4 

.2 

l* 

2.:? 

3. *2 

3.21 

.4 

-  .9 

IP 

2.2 

2.4 

5.5 

2C.4 

55.7 

.2 

l.v* 

:.2i 

Z.li 

•  6 

•  9 

1.9 

2.2 

2.3 

5.8 

iai 

44. 2 

.* 

• . 

2. .2 

3. *2 

2.11 

1.1 

,9 

3.3 

3.3 

3.4 

7.2 

23.2 

51.4 

.* 

21 

2.21 

3.15 

?.»» 

1.9  — 

1.2 

*.* 

4.3 

8.3 

27.3 

*7.0 

•  4 

2  2 

2.21 

».;• 

2.34 

l.€ 

1.9 

9.3 

s.r 

4.3 

M 

23.3 

*5.4 

.* 

25 

2.21 

*.i : 

:.?» 

1  .9 

1.3 

•  .7 

».3 

2.3 

7.4 

2*.* 

*7.2 

1.2 

2* 

2»2r 

>.U 

1.2  - 

1.3 

3.a 

3.3 

2. 2 

4.9 

23.2 

53.* 

2.2 

;« 

2.-1 

•  •24 

2. *2 

t.r 

•  3 

3.2 

:.* 

2.4 

9.2 

23.7 

52. C 

*.3 

Is 

..2* 

i.C 

.7 

3.3 

2.3 

2.3 

4.1 

27.2 

.1.7 

».? 

2’ 

2.i  1 

2. *5 

1.0 

.•  - 

3.4 

2.4 

2.  » 

•  a 

25.2 

•3.3 

12.2 

44 

r 

c 

a 

2 

A 

3 

ft 

5 

6 

7 

• 

, 

:.?i 

*  .  - 
•  a  •* 

•  J4 

9?  .9 

.  w 

t- 

.  • 

.0 

.  - 

2 

1.22 

A**  • 

.  .  ■ 

2’.‘. 

:.i 

•  * 

. : 

. 

•  * 

.  A 

» 

2.33 

2.2  • 

JT.  t 

;•».» 

5t.’ 

#  » 

,  • 

•  T 

.0 

% 

2.22 

•  »*  * 

li.2 

*H.! 

.  w 

•  A 

.  V 

5 

2.33 

’Pft 

:.i  r 

• .  • 

22.3 

Ift.ft 

5.*.: 

:.s 

.  - 

.  • 

% 

a;? 

:.C4 

2.22 

•  *> 

?  .6 

14.7 

j*.# 

2ft.  ft 

.! 

.  2 

.  A 

;.»v 

*.!: 

2.42 

2.3 

11.5 

;b.? 

2’.  ft 

.  • 

* 

a. 2 

T.N 

:pi 

2.4 

2.2 

s  .A 

*.? 

i?.- 

64.2 

.9 

.2 

9 

2.44 

*.*: 

2.74 

2.2 

•  •3 

>.2 

2’.7 

:•  .f 

i!P 

•.ft 

.  A 

1C 

j.i : 

tfl; 

2.5- 

:.t 

2.6 

3.  7 

22.3 

ft*. 7 

2«.5 

. 

a 

2.24 

a?: 

1.3 

».2 

3.2 

!P 

7.4 

!«.( 

rm 

12 

A  .ft  7 

:.* 

.  • 

.! 

A 

.  » 

•  • 

.3 

p 

12.4 

43. C 

12 

2.37 

T  .  *1  - 

2.32 

A 

.* 

.5 

.  ’ 

12.4 

t,2. 6 

!• 

tp: 

’.*7 

•n 

1  .  * 

?.  ft 

2.2 

5.7 

51.* 

P.? 

23 

?.S2 

!.»* 

.  - 

P 

•  .? 

1.7 

•.ft 

27.5 

34.1 

16 

2.23 

2.42 

!.*7 

.  .• 

.7. 

a? 

1 .  * 

1.3 

!?.« 

ft2»2 

.7 

3. .3 

2.iZ 

3.*2 

•  V 

«  1 

1.6 

1.‘ 

4.2 

!».» 

e  *  .  a 

13 

2.4* 

2.45 

3.M 

.2 

•  ‘4 

.? 

#  9 

1  «• 

?.« 

|9.7 

•7.= 

i* 

’.n 

.  •  • 

• .  - 

.* 

01 

1.7 

!  .  * 

2  .ft 

5.7 

52.6 

*•.6 

2m 

:.t* 

a -a 

3.32 

.6 

#  1 

1  •  ft 

!.* 

2.6 

5.7 

25. • 

!:.« 

21 

2. .2 

l.'i 

P 

»  .  “ 

?.  4 

?P 

ir.5 

!ft.R 

i  Z 

2. >3 

2.72 

}.:■> 

A 

. 

i.: 

1.7 

2.1 

2 

4.  ’ 

12. • 

52. * 

2! 

2. *4 

*. 

r.?i 

.7 

i.: 

?  .  * 

?  .ft 

«p 

!'.« 

£•.7 

2» 

2.3* 

r.r-. 

.* 

i.i 

!.“ 

•  . 

f.ft 

21.7 

Sft.6 

25 

a. : 

3.*. 

*••  • 

: . 

!•' 

a: 

2. '* 

2.2 

*.l 

45. ? 

24 

:.~w 

’. a 

3.:. 

•  • » 

3.3 

: .  • 

A  #  A 

4.' 

7 

23 

*.•  4 

2. <*7 

i  •• 

?.2 

:.! 

7.* 

;*.i 

•«.! 

:? 

2.23 

r.  :i 

A  . 

•  .<• 

7.9 

2*  .2 

2' 

i..? 

7. 

i  •  • 

:.r 

’.2 

?.T 

2  #  . 

7.2 

32.2 

pp 

1* 

2.99 

*.7* 

*..*» 

.  .  * 

«  .  ** 

ft  •  • 

:«.4 

17.2 

11 

:.p 

’.’i 

:.? 

•  A 
-  • 

•  • 

5.T 

!■  .a 

S2.« 

22 

3 

4.  •• 

7  .  ft 

*.2 

A  T 

‘P 

:«a 

*1.4 

12 

•  . 

2.4  . 

?.  * 

7. ft 

*•* 

J  .  2 

7.1 

if.* 

.*.7 

24 

i.wC 

••1* 

«. 

1  *  •  4 

1‘lJ 

4  .ft 

;«.? 

13.4 

**-. 

1,* 

A 

*  » 

4 .7 

22.3 

•  1.1 

9 


•C 
•  . 
•  • 
•c 


!*! 


•  9 

•  % 
.• 


. 


TABLE  5  (continued) 
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determined  by  the  amplituae  of  tne  number  in  that  position.  All  the 
graphs  in  the  figure  are  for  1  Hertz.  We  display  the  K..  for  the 

'  J 

moael  S3  in  Table  4  ana  snow  the  E..  at  R  •  244  km  for  this  model 

*  J 

and  the  moael  we  call  S2,  which  has  Q  =  20  to  a  aepth  of  2.8  km.  We 
also  show  tne  E-,  at  1200  km  for  S3.  The  first  and  last  plots 

J 

correspond  to  tabular  data  given  in  Table  5. 

Some  interesting  features  of  the  K..  data  in  tne  table  and 

■  J 

figure  are: 

•  Only  14  modes  contribute  to  1  Hz. 

For  2  Hz,  27  modes  contribute  and  at 
least  35  modes  (the  numDer  computed) 
may  contribute  to  3  Hz. 

•  The  energy  in  tne  funaamental  moae  is 
almost  entirely  in  the  top  600  meters 
of  tne  moael. 

•  The  energy  in  the  higher  modes  is 

almost  entirely  in  the  @  =  3.59 

km/ sec  layers  at  the  base  of  the 
crust. 

•  Energy  in  the  intermediate  crust  is 
restricted  to  low  modes.  At  1  Hz  it 
is  modes  2  ana  3  while  at  3  Hz,  it  is 
modes  2  to  10. 

The  E-  table  is  the  important  one  for  understanding  the 

*  J 

energy  paths  for  the  synthetic  seismograms.  This  quantity  depends 
on  the  Q  model  and  the  range.  Our  observations  about  the  E^  in 
the  taole  and  figure  include  the  following: 

§  Only  20  modes  were  used  to  construct 
the  synthetics.  This  is  enough  to 
compute  1  Hz  and,  apparently,  3  Hz  at 
1200  km.  However,  a  significant 
amount  of  2  Hz  energy  (and  probably 
other  frequencies)  is  lost  by 
truncating  after  20  modes. 

•  The  Q  has  a  dramatic  effect.  At  1  Hz 
the  dominant  energy  is  in  the  5th 
mode  and  is  trapped  at  the  base  of 
the  crust.  The  group  velocity  of 
this  mode  is  3.53  km/sec. 


•  At  2  Hz,  the  modes  9  to  20  dominate 
with  group  velocities  of  3.08  to  3.56 
km/ sec.  At  3  Hz  it  is  modes  12  to  18. 

•  The  effect  of  lowering  the  Q  in  the 

top  2.8  km  (compare  models  S2  and  S3) 
is  to  sharply  attenuate  the  energy  in 
the  low  number  modes  which  have  low 
group  velocities.  This  is  what  we 

saw  in  Figures  3  and  5. 

•  Nearly  all  the  high  frequency  energy 

is  in  the  layer  at  the  base  of  the 

crust.  This  is  easily  seen  in  the 

tabulated  values  of  Ti  as  well  as 

in  the  figure. 

3.2.8  Conclusions 

These  data  ana  the  comparison  of  synthetic  and  observed 

seismograms  suggest  how  the  model  needs  to  be  changed  to  more 

closely  agree  with  the  observations.  The  low  Q  and  gradient  at  the 
surface  were  introduced  to  damp  the  unrealistic  long  period 
fundamental  mode  ringing  that  occurs  with  simple  crustal  models. 

Apparently,  only  a  tnin  low  Q  layer  (thought  of  as  a  laterally 

heterogeneous  weathered  layer)  is  needed. 

For  Lg  propagation,  the  base  of  the  crust  is  most  important. 
A  thick  layer  at  the  base  is  likely  to  give  a  sharp  Lg  onset  at  the 
shear  velocity  of  this  layer,  as  we  see  in  Figure  11.  A  more 

emergent  Lg,  as  seems  to  be  seen  in  the  SALMON  data,  probably 

requires  the  presence  of  some  gradients  in  the  mid-  and  lower- 

crust.  Late  arriving  high  frequency  waves,  say  at  2.9  Km/sec  or 
slower,  can  oe  associated  with  low  number  modes  trapped  near  the 

surface  where  the  shear  velocities  are  low.  However,  high  number 

modes  with  energy  at  the  base  of  the  crust  can  also  have  low  group 

velocities.  In  the  next  section  we  describe  a  new  model  that  gives 

better  agreement  with  the  aata. 

3.3  Lg  FOK  AN  EASTEKN  UNITED  STATES  CKUSTAL  MODEL 

The  parametric  variations  described  in  the  previous  section 

allow  us  to  construct  a  mooel  that  does  give  synthetic  Lg  waveforms 


that  are  qualitatively  like  the  data.  The  new  EUS  model  (called  SI) 
is  listed  in  Taole  6.  It  is  compared  with  the  simple  EUS  crust 
(Table  2)  with  which  we  began  in  Figure  13.  we  see  that  the  velo¬ 
city  model  is  much  the  same  below  4.1  km,  except  that  a  gradient  has 
been  introauced  in  the  mid-crust.  The  Q  is  low  near  the  surface, 

P 

but  increases  to  fairly  high  values  for  depths  below  6.2  km. 

The  phase  and  group  velocity  dispersion  for  the  first  35  modes 
of  model  SI  are  shown  in  Figure  14.  Comparing  to  the  dispersion  for 
Model  S3  in  Figure  7,  we  observe  several  interesting  features 
including: 

•  There  are  no  longer  modes  with  group 
velocities  below  2  km/sec.  This  is 
due  to  the  absence  of  the  very  low 
velocity  surface  layer. 

•  The  higher  mode  Airy  phases  nearly 
all  occur  in  the  2.75  to  3.5  km/sec 
group  velocity  window. 

•  The  band  where  the  group  velocity 
curves  form  a  dense  pattern  is 
somewhat  narrower  than  it  was  for  i3. 

•  There  is  an  indication  of  a  weaker 
set  of  Airy  phases  at  about  2.5 
km/sec. 

Seismograms  were  computed  for  this  model  for  specific 
comparison  to  the  SALMON  data.  The  source  function  of  Mueller  and 
Murphy  (1971)  was  used  with  tne  parameters  for  salt  given  by  Murphy 
(1969).  This  source  function  provides  a  fairly  good  fit  to  the 
measured  near-field  ground  motion  data  when  the  SALMON  yield  (5.3 
KT)  and  depth  (0.83  km)  are  used.  The  data  were  recorded  by  LkSN 
short  period  seismometers  and  the  appropriate  filter  is  applied  to 
the  synthetics. 

Synthetic  seismograms  are  compared  to  the  ooservations  at  two 
ranges  in  Figure  15.  In  the  group  velocity  window  from  3.5  to  2.9 
kmysec,  tne  synthetics  are  like  the  observations  in  a  qualitative 
sense,  particularly  at  the  larger  range.  Tnus,  we  are  encouraged  to 
tnink  that  we  have  done  a  reasonable  job  of  synthesizing  Lg.  Our 


TABLE  6 

EUS  CRUSTAL  MODEL  -  SI 


Depth 

Thickness 

a 

B 

0 

q. 

0.6 

0.6 

3.7 

2.16 

2.10 

20 

2.6 

2.0 

4.55 

2.54 

2.20 

50 

4.1 

1.5 

5.60 

3.14 

2.65 

250 

6.2 

2.1 

6.1 

3.30 

2.85 

400 

13.2 

7.0 

6.3 

3.41 

2.94 

1200 

19.0 

5.8 

6.4 

3.46 

3.00 

1500 

34.0 

15.0 

6.6 

3.59 

3.05 

2000 

00 

GO 

8.1 

4.52 

3.35 
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source  function  are  included. 


15.  (continued) 


The  K..  and  E..  for  the  model  SI  are  plotted  in  Figure  16, 

'  J  'J 

while  the  T.  are  tabulated  in  Table  7.  The  E^  and  I\  were 
computed  for  both  244  km  and  1252  km.  From  the  examination  of  the 
data  in  the  figure  and  table,  we  deduce  the  following: 

•  The  truncation  at  35  modes  appears  to 
be  reasonable.  Most  of  the  energy  is 
in  modes  that  are  near  the  middle  of 
the  band  computed. 

•  At  1252  km  the  high  frequency  energy 
is  nearly  all  in  the  high  Q  layers  at 
the  bottom  of  the  crust.  That  is, 
between  6.2  km  and  34  km  where  the 
shear  velocities  are  3.41  to  3.59 
km/sec. 

•  At  244  km  the  energy  is  more  evenly 
distributed  throughout  the  crust. 

•  The  energy  is  much  more  evenly 
distributed  over  modes  at  closer 
ranges.  At  high  frequencies  and 
larger  ranges,  very  few  modes 
contribute. 

The  tables  indicate  that  the  low  (<  2.8  km/sec)  group 

velocities  are  associated  with  low  number  modes,  which  is  expected 

from  the  plot  in  Figure  14.  From  the  K.  .  plots  we  see  that  the 

J 

energy  in  these  modes  is  trapped  near  the  surface.  Since  the  Q  is 
quite  low  in  the  near  surface  layers  of  our  model,  this  energy  does 
not  propagate.  This  is  seen  in  the  E ■  ■  plots  and  T  table. 

I J  J 

Another  way  to  compare  the  observed  and  synthetic  Lg  phases  is 
in  the  spectral  domain.  In  Figure  17,  we  plot  the  spectrum  of  the 
group  velocity  window  from  3.6  to  2.9  km/sec  for  the  four  vertical 
component  SALMON  observations  in  Figure  12  and  for  four  other 

stations.  The  spectra  from  the  larger  range  (>  1000  km)  stations 

_2 

generally  peak  near  1  Hz  and  fall  off  at  a  rate  between  w  and 

_3 

w  .  At  the  two  closer  stations,  the  spectral  peak  occurs  at 
higher  frequency  and  the  rate  of  falloff  is  about  w  .  The 
signal/noise  is  generally  quite  good  in  the  frequency  band  from  0.7 
to  2.0  Hz.  Outside  this  band  the  signal /noise  is  often  too  small  to 
allow  much  confidence  in  conclusions  about  the  spectral  behavior. 
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synthetics  have  very  little  energy  arriving  before  3.5  km/sec.  As 
long  as  the  mantle  is  represented  by  a  halfspace  (e  -  4.52  km/sec), 
this  will  be  the  case.  Adding  some  layering  in  the  mantle  would 
allow  computation  of  higher  velocity  energy. 

A  more  serious  shortcoming  of  the  model  is  that  there  is  too 
little  late-arriving  high  frequency  energy.  This  is  especially  true 
at  244  Km.  How  can  this  energy  be  included?  The  group  velocity 
curves  in  Figure  14  show  modes  with  stationary  group  velocities 
between  2.6  and  3.0  km/sec.  The  kinetic  energy  distribution  can  be 
used  to  determine  the  portion  of  the  model  that  controls  these  modes. 

The  K..  and  £..  for  the  model  SI  are  plotted  in  Figure  16, 

'  J  '  v 

while  the  are  tabulated  in  Table  7.  The  and  were 

computed  for  both  244  km  and  1252  km.  From  the  examination  of  the 
data  in  the  figure  and  table,  we  deduce  the  following: 

•  The  truncation  at  35  modes  appears  to 
be  reasonable.  Most  of  the  energy  is 
in  modes  that  are  near  the  middle  of 
the  band  computed. 

•  At  1252  km  the  high  frequency  energy 
is  nearly  all  in  the  high  Q  layers  at 
the  bottom  of  the  crust.  That  is, 
between  6.2  km  and  34  km  where  the 
shear  velocities  are  3.41  to  3.59 
km/sec. 

•  At  244  km  the  energy  is  more  evenly 
distributed  throughout  the  crust. 

•  The  energy  is  much  more  evenly 

distributed  over  modes  at  closer 
ranges.  At  high  frequencies  and 
larger  ranges,  very  few  modes 
contribute. 

The  tables  indicate  that  the  low  (<  2.8  km/sec)  group 
velocities  are  associated  with  low  number  modes,  which  is  expected 
from  tne  plot  in  Figure  14.  From  the  K .  .  plots  we  see  that  the 

*  J 

energy  in  these  modes  is  trapped  near  the  surface.  Since  the  Q  is 
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Figure  16b.  The  .  for  structure  SI  are  displayed  for  three  frequencies  at  two  ranges. 

The  maximum  E-jj  is  scaled  to  have  the  value  100.  The  plots  are  then  made 
with  the  shading  indicated  in  the  legend. 


Figure  16b.  (continued) 


TABLE  7 


WEIGHTED  KINETIC  ENERGY  DISTRIBUTION  (Ti) 
BY  LAYER  FOR  THE  MODEL  SI 


I 

H 

BE 

Q 

1  Hz 
Energy 

2  Hz 
Energy 

3  Hz 
Energy 

1 

.6 

2.16 

20 

15.33 

8.49 

6.04 

2 

2.6 

2.54 

50 

21.15 

11.12 

8.58 

3 

4.1 

3.14 

250 

11.92 

9.86 

7.08 

4 

6.2 

3.30 

400 

21.64 

22.22 

11.98 

5 

13.2 

3.41 

1200 

73.20 

79.01 

59.62 

6 

19.0 

3.46 

1500 

55.09 

60.96 

41.86 

7 

34.0 

3.59 

2000 

100.00 

100.00 

100.00 

8 

oo 

4.52 

2000 

1.20 

4.09 

0.09 

1252 

I 

H 

BE 

Q 

1  Hz 
Energy 

2  Hz 
Energy 

3  Hz 
Energy 

1 

0.6 

2.16 

20 

11.45 

3.94 

1.77 

2 

2.6 

2.54 

50 

14.18 

5.82 

2.34 

3 

4.1 

3.14 

250 

7.45 

7.02 

2.62 

4 

65.2 

3.30 

400 

13.04 

14.81 

4.21 

5 

13.2 

3.41 

1200 

54.51 

100.00 

36.57 

6 

19.0 

3.46 

1500 

47.72 

70.62 

13.28 

7 

34.0 

3.59 

2000 

100.00 

98.79 

100.00 

8 

Oo 

4.52 

2000 

1.04 

1.06 

0.03 
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quite  low  in  the  near  surface  layers  of  our  model,  tnis  energy  does 
not  propagate.  This  is  seen  in  the  E-.  plots  and  T.  taole. 

I J  J 

Another  way  to  compare  the  observed  and  synthetic  Lg  phases  is 

in  the  spectral  domain.  In  Figure  17,  we  plot  the  spectrum  of  tne 

group  velocity  window  from  3.6  to  2.9  km/sec  for  the  four  vertical 

component  SALMON  observations  in  Figure  12  and  for  four  other 

stations.  The  spectra  from  the  larger  range  (>  lOOu  km)  stations 

_2 

generally  peak  near  1  Hz  and  fall  off  at  a  rate  between  a>  ana 

_3 

u >  .  At  the  two  closer  stations,  the  spectral  peak  occurs  at 

-3 

higher  frequency  and  the  rate  of  falloff  is  aoout  w  .  The 
signal/noise  is  generally  quite  good  in  the  frequency  Land  from  0.7 
to  2.0  Hz.  Outside  this  band  the  signal/noise  is  often  too  small  to 
allow  much  confidence  in  conclusions  about  the  spectral  behavior. 

In  Figure  18  we  compare  the  synthetic  and  observed  Lg  spectra 
(windowed  between  3.6  and  2.9  km/sec)  from  the  seismograms  plotted 
in  Figure  15.  It  is  not  easy  to  draw  conclusions  from  these 
comparisons,  we  ao  see  that  at  244  km  the  synthetic  spectra  fall 
off  much  more  rapidly  than  the  ooservea  for  frequencies  above  2.6  Hz 
or  so.  The  same  is  true  at  the  larger  distance,  but  the 
observations  are  contaminated  with  noise. 

In  summary,  the  comparisons  of  Figures  15  and  18  lead  to  the 
following  conclusions: 

•  The  Lg  time  signal  is  qualitatively 
similar  to  the  observed  Lg.  The 
onset  of  synthetic  Lg  may  be  a  bit 
more  abrupt  than  is  observed. 

•  Energy  arriving  with  velocities 
greater  than  about  3.6  km/sec  is  not 
modeled. 

t  The  moael  fails  to  account  for  the 

late  arriving  (slower  than  2.9 
km/sec)  high  frequency  energy.  This 
is  associated  with  low  order  modes 
trapped  in  the  top  few  kilometers 
where  Q  is  low. 
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Figure  17.  The  Lg  spectra  (3.6  to  2.9  km/sec)  are  shown  for 
(continued)  observations  of  SALMON  at  eight  LRSM  stations. 

The  spectral  amplitudes  have  been  normalized  and 
lines  of  slope  to-3  are  drawn  on  each  plot.  On  each 
of  the  observed  spectra,  we  plot  a  line  drawn  through 
the  center  of  the  spectrum  of  a  noise  sample  taken 
before  the  Pn  arrival  time.  The  amplitude  of  the 
oscillations  about  this  line  is  similar  to  the  osclll 
tions  in  the  spectra  from  stations  with  R  >  1000  km. 
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Figure  18,  The  Lg  spectra  (within  the  3,6  to  2,9  fcm/sec  group 
velocity  window)  for  the  synthetic  seismograms  In 
Figure  15  are  compared  to  observed  spectra  at  the 
same  range.  The  observed  spectra  are  truncated  by 
a  cosine  filter  above  4  Hz, 


3.4  COMPARISON  OF  OBSERVED  AMO  THEORETICAL  AMPLITUDES 


Thus  far  we  nave  been  entirely  concerned  with  the  waveform  and 
spectral  shape  of  the  Lg  synthetic  seismograms.  A  satisfactory 
model  must  also  give  Lg  phases  that  are  in  reasonable  agreement  with 
observed  amplitude  data. 

There  have  been  numerous  studies  of  the  attenuation  of  the  1 
Hz  Lg  phase  (vertical  component)  in  the  eastern  united  States  (e.g., 
Nuttli,  1S73;  Bollinger,  1973,  1979;  Street,  1976;  uones,  et  al . , 
1977).  These  studies  employ  a  common  method  to  describe  the 
attenuation  rate,  which  was  suggested  by  Nuttli  (1973).  The 
amplitude  is  assumed  to  attenuate  with  distance  as  an  Airy  phase. 
That  is. 


A  *  a-1^3  (sina)-1^  exp(-yA),  (1) 

where  ^  is  some  constant  and  a  is  the  range.  The  attenuation  is 
then  described  by  the  y  which  gives  the  best-fitting  curve  through 
the  observations. 

Nuttli  (1973)  gives  y  *  0.07  deg-*  as  the  value  that  Dest 
characterizes  the  attenuation  of  1  second  vertical  component  Lg  in 
North  America,  east  of  the  Kocky  Mountains.  Bollinger  (1979)  found 
the  0.07  deg-'*'  attenuation  rate  to  be  appropriate  for  epicentral 
distances  from  100  to  700  km  in  the  southeastern  Unitea  States.  At 
longer  distances,  the  rate  is  somewhat  greater  (U.10  deg-1)  for 
some  events.  A  slightly  greater  attenuation  rate  (0.11  deg-*)  was 
found  by  street  (1976)  for  ranges  of  4’  to  30°  in  the  northeastern 
United  States. 

In  Figure  19  we  plot  the  distance  attenuation  curves  for 
several  y  together  with  Lg  amplitudes  measured  from  recordings  of 
sALMun  at  six  stations.  Tnese  recordings  are  displayed  in  Figure 
20.  To  determine  the  Lg  amplitude  measurement,  our  convention  was 
to  measure  the  largest  sustained  amplitude  (several  cycles)  within 
lu  seconds  of  the  group  arrival  time  associated  with  3.5  km/sec 
group  velocity. 


RANGE  (degrees) 


Figure  19,  The  Lg  amplitudes  from  IRSM  recordings  of  SALMON  are 
plotted  versus  range.  Curves  are  plotted  for  several 
apparent  attenuation  rates. 
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Figure  20.  Observed  short  period  recordings  of  SALMON  ai 
minimum  amplitudes  are  indicated  along  with  i 
origin  time  was  16:00).  The  times  for  sever^ 
record. 
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f  SALMON  are  shown  from  six  LRSM  stations.  The  maximum  and 
Tong  with  the  time  of  zero  on  the  time  scale  (the  event 
for  several  apparent  group  velocities  are  indicated  on  each 
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The  attenuation  curves  in  Figure  19  are  forced  to  pass  through 
the  Lg  amplitude  at  EUAL.  The  SALMON  data  appear  to  be  consistent 
with  an  attenuation  rate  ot  about  0.10  to  0.15  deg'1.  However, 
the  data  points  are  restricted  to  two  narrow  distance  bands  (1.2  to 
2.2  degrees  and  9.5  to  12.9  degrees)  so  we  cannot  fit  a  y  curve  very 
confidently.  Further,  the  Lg  attenuation  for  several  of  the  south¬ 
eastern  United  States  events  studied  by  Bollinger  (1979)  changed 
from  0.07  deg'1  for  ranges  less  than  6  degrees  to  0.10°  deg'1 
for  greater  ranges.  Since  the  SALMON  data  are  in  clusters  on  either 
side  of  this  division,  we  are  again  cautioned  about  the  reliability 
of  a  y  determination  from  them. 

In  Figure  21  we  again  plot  the  SALMON  Lg  amplitudes  along  with 
analogous  amplitudes  measured  on  synthetic  seismograms  for  our  model 
SI.  These  synthetics  are  plotted  in  Figure  22  where  we  indicate  the 
Lg  amplitude  measurement. 

The  obvious  conclusion  from  Figure  21  is  that  the  SI  model 
predicts  greater  Lg  attenuation  than  is  observed  in  the  eastern 
united  States.  The  synthetic  data  are  consistent  with  a  Y  or  about 
0.34  deg'1.  The  Y  is  related  to  the  effective  Q  for  the  crust 

P 

by 

Q  *  U7  »  (2) 

where  f  is  the  frequency  and  u  is  group  velocity.  For  1  Hz  energy 
at  3.5  km/sec  group  velocity,  the  effective  Q  is  294. 

The  kinetic  energy  distribution  for  SI  in  Figure  16  and  Table 
7  indicates  that  most  of  the  energy  in  the  synthetic  seismograms  is 
in  the  layers  in  the  lower  crust  where  Q  is  1200  to  2000.  However, 
there  appears  to  be  enough  energy  in  the  low  Q  upper  layers  to  sub¬ 
stantially  reduce  the  effective  Q.  The  effective  Q  can  be  estimated 
by  computing  the  mean  Q"1  weighted  by  the  r.  in  Table  7.  That 


AMPLITUDE  (nm) 


RANGE  (degrees) 

Figure  21.  The  synthetic  Lg  amplitudes  from  the  SI  synthetic  seismo¬ 
grams  are  plotted  with  the  SALMON  data. 
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Figure  22.  (continued) 


where 


(3) 


n 


i«l 


This  effective  Q  should  not  be  confused  with  the  Rayleigh  Q  of  a 
particular  mode,  but  it  is  an  estimate  of  the  total  Q  for  all  modes 
which  contribute  to  the  seismograms  at  a  particular  frequency.  The 
Qgff  values  computed  from  the  for  SI  in  Table  7  are  listed 
below: 


R  -  244  km 

R  *  1252  km 

1  Hz 

208 

238 

2  Hz 

328 

540 

3  Hz 

356 

655 

Tne  effective  Q  associated  with  the  y  fit  in  Figure  21  is  294,  which 
falls  oetween  the  1  and  2  Hertz  values  in  the  above  listing.  The 
dominant  frequency  where  the  Lg  amplitudes  are  measured  in  tne 
seismograms  in  Figure  22  is  in  this  range. 

The  good  agreement  between  the  effective  Q  uetermined  oy 
huttli's  ( 1973 J  method  in  Figure  21  and  the  Qgtf  uetermined  from 
the  mean  Q-1  weighted  by  the  amount  of  kinetic  energy  in  each 
layer  is  quite  an  interesting  result.  Note  that  the  Qe^ 
increases  with  range.  This  means  that  the  path  for  the  dominant 
energy  is  deeper  for  larger  ranges,  which  is  in  accordance  with  our 
expectation.  The  also  increases  with  frequency.  This  is 

because  the  higher  frequency  energy  at  shallow  depths  is  severaly 
attenuated  and  contributes  very  little  energy  to  the  seismogram. 

The  increase  in  with  range  suggests  that  the  y  should 

decrease  with  range.  In  fact,  the  amplitude  data  in  Figure  21  show 
this  trend.  The  two  smallest  range  plots  are  fit  with  y  *  u.42, 
which  implies  a  1  Hz  effective  Q  of  238,  while  the  two  larger  range 
points  are  fit  with  y  «  0.32,  which  implies  an  effective  Q  of  312. 


3.5  PARAMETRIC  VARIATIONS  OF  THE  MODEL  SI 
3.5.1  Introfluction 

The  model  SI  leads  to  synthetic  Lg  that  resemble  observed  Lg 
in  the  group  velocity  window  from  3.5  to  2.9  km/sec.  The  associated 
spectra  are  also  in  reasonable  qualitative  agreement  with  the  data, 
but  there  are  several  data  features  that  are  not  modeled  very  well 
by  the  synthetic  seismograms.  These  are: 

The  synthetics  include  too  little  late  arriving 
high  frequency  energy. 

The  synthetic  Lg  spectra  fall  off  too  rapidly 
at  high  frequencies,  especially  at  the  closer 
range. 

The  synthetic  Lg  amplitude  attenuation  is  much 
greater  tnan  observed  in  the  eastern  United 
States. 

In  this  section  we  show  the  results  from  several  mooels  in 
wnicn  we  attempted  to  improve  the  synthetic  seismogram  match  to  the 
data  characteristics  mentioned  aoove.  Most  of  our  attention  is  on 
changing  the  Q  model  to  reduce  the  attenuation,  since  the  listed 
model  deficiencies  all  seem  to  indicate  that  bl  has  too  much  attenu¬ 
ation.  However,  we  are  constrained  by  the  requirement  for  low  Q 
near  the  surface  to  reduce  the  amplitude  of  the  fundamental  mode,  as 
was  discussed  in  Section  3.2.  Also,  we  saw  in  Figure  21  that  the 
agreement  of  observed  and  synthetic  Lg  amplitudes  is  really  not  too 
bad  as  far  as  absolute  values  are  concerned.  Simply  raising  the  Q 
to  match  the  observed  amplitude-distance  attenuation  will  not  be 
satisfactory  because  it  will  give  synthetic  seismogram  amplitudes 
that  are  huge  compared  to  the  observed  SALMON  amplitudes. 

In  this  section  we  will  first  describe  more  carefully  the 
source  functions  and  instrument  response  used  in  in  constructing  the 
synthetic  seismograms,  while  there  may  be  some  errors  in  these 
factors,  especially  the  source  specification,  they  are  not  likely  to 
influence  the  gross  features  we  are  looking  for  in  our  comparison 
witn  data. 


65 


FREQUENCY  (HZ) 


Figure  23.  Spectre!  amplitude  of  the  Instrument  response  and 
source  function  used  to  construct  the  synthetic 
seismograms.  There  appears  to  be  an  anomalous 
point  In  the  Instrument  response  near  1,6  Hz, 
but  this  should  have  little  effect. 
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FILE=ELASTIC  SI 

STATION:  MOD=35/EVENT :  8=.83,R=244 
TSTART=  0.  :  0.:  -29.0000 


Figure  24.  Synthetic  seismograms  for  the  crustal  model  SI  with  Infinite  Q  are  compared 
to  observations  of  SALMON.  Zero  time  on  the  synthetics  Is  correctly  aligned 
with  the  observations. 


VOTO  16.  3.  30. 


TABLE  8 


WEIGHTED  KINETIC  ENERGY  DISTRIBUTION  (I^)  BY 
LAYER  FOR  SI  WITH  INFINITE  Q 


H 

BE 

Q 

1  Hz 
Energy 

2  Hz 
Energy 

3  Hz 
Energy 

0.60 

2.16 

00 

100.00 

100.00 

2.60 

2.65 

00 

61.50 

18.21 

12.92 

3.14 

00 

7.20 

3.71 

2.59 

6.20 

3.30 

00 

4.06 

1.11 

1.52 

13.20 

3.41 

00 

9.44 

1.35 

2.68 

19.0 

3.46 

00 

6.64 

1.00 

2.02 

34.00 

3.59 

oo 

11.84 

2.13 

1.37 

35.00 

4.52 

00 

0.15 

0.07 

0.00 
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SPECTRAL  AMPLITUDE 


Figure  26.  The  Ig  spectra  are  plotted  at  two  ranges  for  the 
model  SI  with  infinite  0.  These  spectra  are  trun¬ 
cated  with  a  cosine  filter  above  4  Hz. 
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Since  tne  Q  model  is  the  main  issue,  we  look  at  the  synthetics 
for  infinite  g  as  a  bound.  We  then  compute  seismograms  for  several 
Q  models  that  might  be  thouyht  to  improve  agreement  with  the  data. 
We  also  look  at  one  model  in  which  the  velocity-depth  profile  is 
somewhat  different  from  SI,  though  the  crustal  thickness  is  held 
fixed.  Our  conclusions  based  on  these  parameter  variations  are 
listed  in  Section  3.6. 

3.5.2  Source  and  Seismometer  Contributions 

In  Figure  23  we  plot  the  spectrum  of  the  source  function  and 
seismometer  response  used  for  the  synthetic  seismogram 
calculations.  The  source  is  presented  in  terms  or  the  spectral 
amplitude  of  the  reduced  velocity  potential.  This  is  the  SALMON 
source  function  proposed  by  Murphy  (1969)  whicn  gives  good  fit  to 
tne  observed  near-field  ground  motions. 

The  source  spectrum  is  nearly  that  of  a  delta  function  over 
the  relevant  frequency  band.  The  actual  source  function  would  have 
to  oe  very  strongly  peaked  at  frequencies  greater  than  2.5  Hz  to 
bring  the  observed  and  theoretical  spectra  into  closer  agreement. 
Tne  near-field  data  (apringer,  et  al . ,  1968)  give  no  indication  that 
this  is  the  case. 

3.5.3  An  Elastic  Model 

What  is  the  effect  of  the  Q  model  we  have  assumed?  In  Figure 
24  we  show  the  seismograms  for  the  structure  Si  with  infinite  Q. 
That  is,  the  model  is  elastic.  Later  arriving  high  frequency  enwgy 
is  present,  but  the  seismogram  is  dominated  by  the  sinusoidal 
fundamental  mode  waveform.  This  can  be  clearly  seen  in  the  E.j 
plots  in  Figure  25  and  r.  listing  in  Table  8.  Nearly  all  the 
energy  at  high  frequencies  is  in  the  fundamental  mode  and  is  trapped 
in  the  top  600  meters.  Even  at  1  Hz  only  the  top  two  layers, 
extending  to  a  depth  of  2.6  km,  are  important. 

The  Lg  spectra  for  the  elastic  model  seismograms  are  plotted 
in  Figure  26.  Comparing  to  the  observed  spectra  in  Figure  17,  we 
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see  that  the  elastic  model  does  lead  to  spectra  with  about  the  right 
high  frequency  falloff,  though,  the  high  frequencies  are  actually 
too  dominant  in  these  spectra.  An  elastic  model  is  by  no  means 
thought  to  be  realistic,  but  does  indicate  that  higher  Q  will  bring 
the  high  frequency  part  of  the  synthetic  spectra  in  closer  agreement 
with  observations. 

3.5.4  Tne  Effect  of  Variations  in  the  Q  Model 

In  Taole  9  we  list  four  Q  moaels  for  which  seismograms  will  be 
compareci.  The  first  is  SI  from  Table  6  ana  all  the  moaels  have  this 
velocity  profile.  The  next  model,  S5,  has  higher  Q  in  the  mid-crust 
and  much  higher  Q  in  the  near-surface  layers.  The  Model  S6  has  a 
very  low  Q  in  the  top  600  meters  ana  is  essentially  elastic  at  the 
base  of  the  crust.  In  the  mid-crust  it  is  a  high  Q  nioael.  An 
extreme  example  of  a  hign  0  moael  is  37,  for  which  only  low  Q  values 
are  in  the  top  2.6  km.  Synthetic  seismograms  are  compared  for  these 
four  models  in  Figure  27. 

The  Lg  amplitudes  of  these  seismograms  are  listed  in  Table  lu 
along  with  the  observed  SALMOH  amplitudes  from  stations  at  tne  same 
range.  Models  o5  and  S7  have  amplitudes  very  close  to  those 
observed  at  the  large  distance,  and  S6  is  not  too  mucn  different. 
However,  tne  nigner  Q  models  all  have  amplitudes  much  larger  than 
those  observed  at  the  244  km  range.  This  is  because  the  attenuation 
parameter  y  remains  too  large. 

To  quantify  the  attenuation  we  again  look  at  tne  energy 
distribution  in  these  models  in  Figure  28  and  Table  li.  The  y  for 
these  moaels  can  be  computed  from  the  in  Table  11  using  (2) 
with  the  Qg1;f  computed  from  (3)  in  Section  3.4.  This  gives  a  y 
for  a  specific  range,  we  know  that  y  decreases  witn  range  because 
the  proportion  of  energy  in  the  higner  Q  deep  layers  is  greater. 

We  compute  tne  y  for  S5,  S6  and  S7  at  R  -  1242  km  and  for  f  » 
1.5  Hz.  The  at  1.5  Hz  is  taken  to  be  the  mean  of  the  Q  ^ 

values  for  1  and  2  Hz,  computed  from  (3)  using  the  r.  in  Table 
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TABLE  9 


Q  MODELS  FOR  THE  STRUCTURE  OF  TABLE  3  MODEL 


Depth 

8 

SI 

S2 

S6 

0.6 

2.16 

20 

50 

20 

2.6 

2.54 

50 

350 

200 

4.1 

3.14 

250 

400 

400 

6.2 

3.30 

400 

500 

500 

13.2 

3.41 

1200 

2000 

3000 

19.0 

3.46 

1500 

2000 

4000 

34.0 

3.59 

2000 

2000 

5000 

00 

4.52 

2000 

2000 

5000 

S7 

20 

200 

2000 

5000 

10000 

10000 

10000 

10000 
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285  SECONDS 


Figure  27.  (continued) 


TABLE  10 


Lg  AMPLITUDES  FOR  THE  SEISMOGRAMS  IN  FIGURE  27 
Lg  Amplitude  (run) 

Model  R  »  244  km  R  -  1252  km 


SI 

5250 

58 

$5 

11140 

266 

S6 

8930 

167 

S7 

11590 

265 

Observed  Amplitude 


distributions  are  plotted  for  the  models  S5  and  S6  at  two  frequencies 
:s  for  S7  are  nearly  the  same  as  those  for  S6. 


WEIGHTED  KINETIC  ENERGY  DISTRIBUTION  (r<)  IN  THREE  Q  MODELS 
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11.  The  y  computed  this  way  are  0.21  deg-*  for  S5  and  S 7  and  0.26 
deg-*  for  $6.  These  should  be  compared  to  r  *  0.34  deg-*  for  SI 
which  we  found  in  Section  3.4.  Our  confidence  that  these  are  the 
correct  y  characterizing  these  four  structures  is  enhanced  by  the 
fact  that  they  almost  perfectly  predict  the  relative  Lg  amplitudes 
at  K  >  1252  km  listed  in  Table  10. 

The  final  comparison  is  that  for  the  Lg  spectra  and  is  shown 
in  figure  29.  These  spectra  are  remarkably  similar.  Since  all  the 
seismograms  were  computed  with  the  same  velocity  model,  ana  the  Lg 
is  mostly  associated  with  propagation  near  the  base  of  the  crust 
where  the  Q  is  high  in  every  model,  this  should  not  be  too 
surprising.  The  lowest  Q  model  is  SI,  and  this  is  evidenced  in  the 
high  frequency  part  of  the  spectrum,  especially  at  1252  km. 

3.5.5  Lg  Synthetic  Seismograms  for  an  Altered  Velocity  Model 

Before  listing  our  conclusions  from  comparing  the  Q  models 
listed  in  Table  9,  we  examine  the  seismograms  from  a  model  with  a 
somewhat  different  velocity-depth  profile.  This  is  the  mooel  S4 
listed  in  Table  12.  The  b  and  Q  for  this  model  are  compared  to 
those  for  SI  in  Figure  30.  We  see  tnat  the  two  models  have  the  same 
crustal  thickness.  The  high  velocity  crustal  granites  are  closer  to 
the  surface  in  S4  ana  a  small  gradient  has  Deen  added  at  trie  base  of 
the  crust.  The  has  higher  Q  in  the  top  4  km  and  lower  Q  in  the 
mid-crust  between  4  and  12  km. 

The  phase  and  group  velocity  curves  for  S4  are  shown  in  Figure 
31.  Comparing  to  those  for  51  in  Figure  12,  we  see  that  there  are 
relatively  few  modes  with  group  velocities  less  than  3  km/sec,  oue 
to  there  being  less  gradient  near  the  surface. 

The  seismograms  for  this  model  are  shown  in  Figure  32,  where 
they  are  compared  to  the  SALMUN  data  and  to  the  seismograms  from 
51.  The  Lg  has  a  more  abrupt  onset  for  S4.  It  is  also  longer 
period.  This  can  be  seen  in  Figure  33  where  we  compare  Lg  spectra 
for  the  two  models.  The  lower  Q  in  the  mid-crust  is  apparently  the 
cause  for  this  frequency  shift.  Also,  as  we  have  seen  before,  the 
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Figure  30.  Comparison  of  models  SI  and  S4 


FREQUENCY  (HZ) 


Figure  31.  The  first  35  modes  for  the  model  S4 


Figure  32.  Thirty-five  synthetic  seismograms  for  the  models  SI  (Figure  15)  and  S4  are  compared 
to  SALMON  observations  from  two  stations.  The  peak  amplitude  In  microns  is  listed 
at  the  left  of  each  synthetic.  The  time  scales  are  aligned  at  3.5  km/sec. 


32.  (continued) 


SPECTRAL  AMPLITUDE 


relatively  hign  Q  near  the  surface  in  $4  admits  a  large  fundamental 
mode  contribution  with  a  sinusoidal  waveform  not  seen  in  the  data. 

In  Figure  34  we  show  the  K..  and  E . .  distributions  for 

’  J  *  J 

this  model.  The  iv  are  tabulated  in  Table  13.  From  the  K •  ^  we 
see  immediately  that  the  added  layer  at  the  base  of  the  crust  will 

not  have  a  strong  effect.  The  e  =  3.59  km/sec  layer  still  dominates 

the  kinetic  energy  distrioution  for  the  nigher  modes. 

The  Lg  amplitudes  for  the  seismograms  in  Figure  32  are  6550  nm 
at  k  *  244  and  92  nm  for  R  =  1252  km.  These  are  between  the 
amplitudes  for  SI  and  S6  in  Table  10.  That  is,  like  those  from  all 
the  other  models,  they  fall  off  much  too  rapidly  between  these  two 
ranges.  The  y  compute-  using  (2)  and  (3)  and  the  r,  at  k  *  1252 

km  in  Table  13  is  0.34  deg  .  This  is  the  1.5  Hz  y  computed  tne 
same  way  we  determined  the  y  for  models  55,  S6  and  :>7.  This  value 
is  quite  close  to  tne  y  we  aerived  for  the  moael  SI  in  Section  3.3. 

3.6  CONCLUSIONS 

In  Section  3.2  we  found  that  some  low  velocity,  low  Q  material 
is  needed  in  the  near  surface  (<  1  km)  region  to  prevent  domination 
by  an  unrealistic  fundamental  mode  harmonic  wavetrain.  As  we  varied 
the  0  moaels  in  the  last  section,  we  saw  that  in  several  cases, 

especially  S4  and  55,  the  near  surface  Q  was  about  as  high  as  it 
could  be  ano  not  have  the  fundamental  mode  become  too  large  to  be 
realistic.  Yet  none  of  the  models  include  tne  high  frequency  energy 
arriving  after  2.9  km/sec  that  is  seen  on  the  observed  SALMON 
recordings.  The  higher  Q  models  do  have  larger  amplitudes  at  lower 
velocities;  for  example,  compare  SI  and  57  at  k  *  1252  km  in  Figure 
27,  but  still  do  not  have  much  energy  after  2.9  km/sec. 

The  spectral  comparisons  are  in  Figures  29  and  33.  Recall 
from  Figure  17  that  the  data  have  a  high  frequency  slope  near  u> 
at  the  closest  stations  and  seem  to  fall  off  at  a  rate  between 
and  u>-J  at  the  larger  ranges  \some  uncertainty  is  introduced  at 
the  larger  ranges  by  the  influence  of  seismic  noise).  The  high 


Figure  34.  The  K. .  and  E.,  distributions  are  plotted  for  the  model  S4 
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frequency  character  of  the  data  is  not  reproauced  by  the  synthetic 
seismograms  for  any  of  the  models. 

The  high  Q  models  $5,  S6  and  $7  were  selected  to  reduce  the  y 
to  a  value  closer  to  that  observed  in  the  eastern  tmiteu  States 
(section  3.4).  The  models  SI  and  S4  have  y  values  of  about  0.34 
deg-^,  much  higher  than  the  0.07  to  0.10  deg-^-  values  observeu. 
The  mouels  S5  and  S7  do  reduce  the  y  to  0.21  deg“^,  which  is  an 
improvement,  but  still  much  too  high.  These  two  models  are  quite 
different,  even  though  they  result  in  the  same  Lg  amplitudes  (Table 
10),  which  match  the  SALMON  observation  at  1252  km,  but  are  a  factor 
of  4  larger  at  244  km.  The  S5  model  has  relatively  high  Q  values  in 
tne  top  2.6  km,  which  allows  too  much  fundamental  mode  energy  at  the 
closer  range  (Figure  27).  The  S7  model  is  elastic  below  13  km  and 
has  very  high  Q  below  4  km.  However,  it  also  has  a  very  low  Q  near 
surface  layer  which  does  succeed  in  damping  the  fundamental  mode. 

Thus,  we  are  left  with  the  deficiencies  we  discussed  in 
connection  with  Moael  SI  ana  Section  3.3.  These  are: 

1.  The  synthetic  seismograms  have  too  little  high 
frequency  energy  at  velocities  less  than  2.9 
km/sec. 

2.  The  synthetic  Lg  spectra  (windoweo  between  3.5 
and  2.9  kin/sec)  have  too  little  high  frequency 
energy. 

3.  Tne  y  for  the  synthetic  seismograms  is  much 
larger  than  observed. 

Our  parametric  studies  are  certainly  far  from  complete.  *Je 
have  not  yet  looked  at  much  different  velocity-depth  profiles,  par¬ 
ticularly  with  a  different  crustal  thickness.  There  are  many  vari¬ 
ations  on  the  Q  model  that  also  could  be  tried.  However,  these 
parametric  studies  are  enough  to  strongly  indicate  that  we  will  not 
be  able  to  remove  the  listed  deficiencies  with  plane-layered,  con¬ 
stant  Q,  laterally  homogeneous  models.  The  late-arriviny  high  fre¬ 
quency  energy  may  be  largely  aue  to  scattering  and  multipathing. 
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3ut  the  most  important  upgrading  of  the  models  is  likely  to  be  the 
introduction  of  a  frequency  dependent  Q.  Recent  results  (Mitchell, 
1980)  indicate  that  this  is  important  to  match  amplitude  attenuation 
data  in  the  eastern  united  States. 


IV.  THE  DEPENDENCE  OF  Lg  ON  SOURCE  DEPTH 


4.1  INTRODUCTION 

An  important  question  is  how  strongly  the  Lg  phase  depends  on 
the  source  depth.  To  address  this  question,  we  compute  synthetic 
seismograms  in  the  model  SI  of  section  3.3  Table  6.  While  this 
model  does  not  give  seismograms  that  match  the  observed 
characteristics  of  Lg  in  all  respects,  it  should  be  good  enough  to 
indicate  the  gross  effects  of  varying  source  depth.  Our  initial  set 
of  seismograms  were  done  with  a  center  of  dilatation  source.  This 
is  not  a  realistic  source  representation  for  depths  below  a  few 
kilometers,  but  does  lead  to  some  interesting  results  that  are  worth 
mentioning.  In  fact,  these  results  show  that  the  center  of 
dilatation  source  gives  an  unrealistic  estimate  for  the  depth 
dependence  of  the  higher  modes  that  control  Lg. 

The  main  calculations  to  determine  the  depth-dependence  of  Lg 
were  done  with  a  point  double-couple  source  at  the  three  basic 
orientations  (strike-slip,  vertical  dip-slip  and  45°  dip-slip)  which 
can  be  combined  to  give  the  solution  for  a  double-couple  with 
arbitrary  orientation  (Burridge,  et  aj_. ,  1964). 

4.2  DEPTH-DEPENDENCE  OF  Lg  FOR  A  CENTER  OF  DILATATION  SOURCE 

Synthetic  seismograms  were  computed  at  seven  depths  from  1  to 

14  kilometers  in  the  crustal  model  SI.  These  seismograms  are  shown 

in  Figure  35.  To  compare  amplitudes,  the  source  was  scaled  to  a 

22 

common  moment  (1.7  x  10  dyne-cm) .  Since  the  moment  of  a  center 
of  dilatation  source  is  given  by  (e.g.,  Hudson,  1969) 

Mq  .  4»pa2  4'00, 

where  f^is  the  static  value  of  the  reduced  displacement  potential, 
we  fix  f,*,  to  be  2960  m^  and  scale  by  pa^.  The  peak-to-peak 
amplitude  is  indicated  at  the  left  of  each  record. 
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Figure  35.  Synthetic  seismograms  are  shown  for  the  SALMON  source  at  seven  focal  depths. 

The  amplitude  scale  at  left  is  in  nanometers.  The  source  moment,  which  is  pro¬ 
portional  to  pa2,  is  held  constant.  The  source  region  a  is  indicated  with  each 
record. 
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The  seismograms  have  a  very  strong,  but  irregular,  dependence 
on  depth.  This  is  not  a  general  characteristic  of  Lg,  but  occurs 
only  for  a  center  of  dilatation  source.  The  excitation  of  the  high 
frequency,  higher  mode  energy  dies  off  rapidly  away  from  layer  inter¬ 
faces  for  this  P  wave  source.  In  the  model  SI,  there  are  interfaces 
at  depths  of  6.2  and  13.2  km.  Thus,  the  8  and  11  km  depths  are  near 
the  center  of  this  layer,  while  the  6  and  14  km  depths  are 
relatively  near  the  interfaces. 

From  the  seismogram  comparison  in  Figure  35  it  is  obvious  that 
the  depth  has  a  much  stronger  effect  on  the  excitation  of  high 
frequency  energy.  This  is  also  seen,  of  course,  in  the  spectra  of 
these  seismograms.  We  plot  these  in  the  form  of  what  we  Cull 
"smoothed  Lg  spectra."  These  are  obtained  by  first  computing  the 
Fourier  spectra  of  the  synthetic  seismograms  in  the  3.6  to  2.9 
km/sec  group  velocity  window.  These  spectra  are  then  smoothed  by 
computing  a  thirty-one  point  moving  average.  A  typical  application 
of  the  smoothing  is  shown  in  Figure  36.  The  smoothed  Lg  spectra  for 
the  seismograms  in  figure  35  are  shown  in  Figure  37. 

4.3  DEPTH-DEPENDENCE  FOR  A  DOUBLE-COUPLE  SOURCE 

To  investigate  the  depth-dependence  of  a  double-couple  source, 

synthetic  seismograms  were  computed  at  various  depths  in  the  model 

SI  of  Section  3.3,  Table  6.  The  source  is  a  point  double-couple 

25 

with  a  fixed  moment  of  10  dyne-cm.  The  dislocation  time  history 
is  a  step  function,  so  the  source  has  no  corner  frequency.  Seismo¬ 
grams  were  computed  for  this  source  oriented  to  represent  strike- 
slip,  vertical  dip-slip  and  45*  dip-slip  faulting.  The  solution  at 
any  other  fault  orientation  can  be  obtained  by  summing  solutions  for 
these  three  basic  orientations  (Burridge,  et  al_.,  1964). 

Synthetic  seismograms  which  demonstrate  the  depth-dependence 
of  Lg  are  shown  in  Figures  38  to  43.  First,  in  Figures  38,  39  and 
40,  we  show  the  seismograms  at  seven  depths  from  1  to  14  km  for  all 
three  source  orientations.  The  orientation  of  the  source  is  given 
by  dip  («)  and  slip  (x)  angles  and  the  source-station  azimuth  is  e. 
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DEPTH  (KM) 


Figure  37.  Smoothed  Lg  spectra  for  the  center  of  dilatation  source  of 
Figure  35.  These  spectra  are  cosine  tapered  to  zero  between 
4  and  5  Hertz. 
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Figure  38.  Synthetic  seismograms  for  a  strike-slip  (X  =  0°,  5  =  90°, 

0  -  30°)  double-couple  with  a  step  dislocation  time  history 
with  moment  10“  dyne-cm.  The  range  is  1252  km  and  the  LRSM 
short  period  instrument  response  is  included.  An  amplitude 
measurement  was  made  at  the  cycle  marked  with  an  arrow  which 
is  within  5  seconds  of  time  associated  with  a  group  velocity 
of  3.5  km/sec. 
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Figure  39.  Synthetic  seismograms  for  a  normal  dip-slip  (X  =  90°,  5  =  90 
0  *  22.5°)  double-couple.  The  range  is  1252  km  and  the  LRSM 
short  period  response  is  included. 
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Figure  40.  Synthetic  seismograms  for  a  45°  dip-slip  (\  *  90°,  6  =  45°, 

9  =  22.5°)  double-couple.  The  range  is  1252  km  and  the  LRSM 
short  period  response  is  included. 
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These  seismograms  have  been  filtered  by  the  LRSto  short  period 
instrument  response.  For  all  three  orientations  we  show  the 
seismograms  at  a  range  of  1252  km.  For  the  strike-slip  orientation 
(Figure  36),  we  also  show  the  seismograms  at  a  range  of  244  km. 

To  demonstrate  the  dependence  over  a  wider  depth  variation,  a 
different  set  of  seismograms  is  shown  in  Figures  41,  42  and  43.  In 
tnis  case,  the  range  is  1000  km  and  the  wWSSN  short  period 
instrument  response  is  included.  The  depths  are  5,  15  and  25  km. 

There  seems  to  be  little  systematic  dependence  of  the 
waveforms  on  depth.  The  duration  of  the  normal  dip-slip  synthetics 
appears  to  decrease  with  depth.  Otherwise,  the  primary  depth  effect 
is  that  the  amplitude  changes.  The  Lg  amplitude  was  measured  on 
each  record  and  these  are  plotted  versus  source  depth  in  Figure  44. 
In  this  case,  we  measured  the  maximum  peak-to-peak  amplitude  witnin 
five  seconds  of  3.5  km/sec.  This  phase  is  marked  with  an  arrow  on 
each  seismogram.  The  amplitudes  plotteu  in  Figure  44  are  the  raw 
amplitudes  from  the  seismograms;  that  is,  no  corrections  were  made 
for  frequency  dependent  instrument  response. 

The  depth  of  tne  layer  interfaces  is  shown  in  the  figure  along 
with  tne  Q.  at  each  depth.  There  is  little,  if  any,  indication 

P 

that  the  amplitude  depends  on  the  layer  in  which  the  source  occurs. 

To  see  tne  depth  dependence  over  tne  full  range  of  depths,  we 
need  to  scale  the  amplitudes  from  the  two  ranges  and  seismometers  so 
they  can  be  plotted  together.  To  do  so  we  scaled  tne  WW5SN,  R  = 
1000  km,  Lg  amplitudes  so  they  fit  on  a  line  drawn  between  the  4  and 
6  km  depth  points  on  the  amplitude  plots  for  the  LfO>M,  R  *  1252  km 
synthetics.  The  resulting  composite  curves  are  shown  in  Figure  45. 
There  is  a  small  break  in  the  curve  between  the  14  and  15  km  depth 
points,  which  is  probably  due  to  the  range  anu  seismometer  differ¬ 
ences,  but  tnis  does  not  disturb  the  main  trends  apparent  in  the 
plot. 

we  conclude  tnat  tne  Lg  amplitude  decreases  with  depth  at 
nearly  the  same  rate  for  the  strike-slip  anu  45*  oip-slip  source 
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Figure  41.  Synthetic  seismograms  for  a  strike-slip  (X  *  0°,  <5  *  90°, 

8  *  30°)  double-couple.  The  range  is  1000  km  and  the  WWSSN 
short  period  response  is  included. 
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Figure  42.  Synthetic  seismograms  for  a  normal  dip-slip  (X  *  90°,  f  =  90°, 
0  =  30°)  double-couple.  The  range  is  1000  km  and  the  WWSSN 
short  period  response  is  included. 
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Figure  43.  Synthetic  seismograms  for  a  45°  dip-slip  (>  =  90°,  6  =  45°, 
6  *  30°)  double-couple.  The  range  is  1000  km  and  the  WWSSN 
short  period  response  is  included. 


:l] 


SctMOMter 


44.  Lg  amplitudes  from  the  synthetic  seismograms  of  Figures  38  to  43  are  plotted  versus 
source  depth.  The  Interfaces  are  Indicated  with  a  dashed  line  and  the  Q  In  each  la 
Is  listed  along  the  right  border  of  the  plot.  The  amplitude  scale  for  the  R  =  244 
synthetics  Is  at  the  top. 


comparable  and  replotted.  The  scaled  points 
rked.  The  format  is  like  Figure  44,  except  the 
eft  border. 


orientations.  This  is  a  factor  of  four  to  five  decrease  in 
amplitude  between  1  and  25  km  depths.  For  vertical  dip-slip 
faulting,  the  Lg  amplitude  increases  with  depth  all  the  way  to  the 
aeepest  (25  km)  computed.  The  rate  of  increase  is  a  factor  of  five 
of  six  over  the  depths  plotted. 


Ill 
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